Mon. Not. R. Astron. Soc. 000,[T]|2^(0000) Printed 5 February 2008 (MN MgX style file v2.2) 



X-ray Scaling Properties of Early-type Galaxies 

Ewan O' Sullivan*, Trevor J. Ponman, Ross S. Collins 

School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT 



Accepted 2003 ??. Received 2002 ??; in original form 2002 ?? 



m 
o 
o 

(N 



> 

m 
in 



o 
cn 
o 

:^ 

Oh; 

6 



X 



ABSTRACT 

We present an analysis of 39 X-ray luminous early-type galaxies observed with the ROSAT 
PSPC. Using multi-component spectral and spatial fits to these data we have measured halo 
abundance, temperature, luminosity and surface brightness profile. We compare these mea- 
surements to similar results from galaxy groups and clusters, fitting a number of relations 
commonly used in the study of these larger objects. In particular, we find that the (t:Tx re- 
lation for our sample is similar to that reported for clusters, consistent with (3spec — 1, and 
that the Lx'Tx relation has a steep slope (gradient 4.8±0.7) comparable with that found for 
galaxy groups. Assuming isothermality, we construct 3-dimensional models of our galaxies, 
allowing us to measure gas entropy. We find no correlation between gas entropy and system 
mass, but do find a trend for low temperature systems to have reduced gas fractions. We con- 
clude that the galaxies in our sample are likely to have developed their halos through galaxy 
winds, influenced by their surrounding environment. 
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1 INTRODUCTION 

Early-type galaxies have been known to possess large halos of hot 
gas since the detection by Einstein of X -ray emi ssion from the el- 
liptical population in the Virgo cluster jiForman et alJll979l) . Suc- 
cessive generations of X-ray observatories have been used to ob- 
serve these galaxies, and the advent of XMM-Newton and Chan- 
dra has allowed the complex nature of their emission to be studied 
in detail. Most of the work in this area has focused on the vari- 
ous sources of emission within early-type galaxies (hot gas. X-ray 
binaries, AGN) and on the surprisingly complicated relation be- 
tween optical and X-ray luminosity. However, at a quite fundamen- 
tal level, early-type galaxies resemble the groups and clusters in 
which they typically reside. Simulations of dark ma tter halos sug- 
gest t hat they have similar profiles at all mass scales iNavarro et alJ 
Il997ft . If we consider clusters, groups and galaxies as potentials 
containing hot gas, we might expect the properties of the halos to 
be similar across a wide range of masses. 

A comparison between halos on different scales becomes 
increasingly interesting considering the importance of entropy 
changes in governing the behaviour of group and cluster halos. 
Observations of galaxy groups have demonstrated that these sys- 
tems do not behave as might be expected from scaling clusters, 
but instead require non-gravitation al processes. One of the clear- 
est signs of this is the entropy floor jPonman et alll999^ . Whereas 
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in more massive systems gas entropy scales with the total mass, 
in groups it appears to reach a roughly constant minimum level. 
A number of models have been put forward to explain this be- 
haviour, including raising the entropy through the injectio n of en- 
ergy by AGN jWu et al l boOO) or star formation I Ponm an et alJ 
Il999|). or th rough the radiative cooling and removal of low entropy 
gas iMuaiiwong et al. 2001.) . It is notable that for any of these pro- 
cesses, the origin of the entropy rise would be related to the galax- 
ies in the system. As early-type galaxies possess their own halos, 
we might expect to see evidence of these processes in their X-ray 
properties, and for the effect to be strongest in these systems, owing 
to their position at the bottom of the mass scale. 

The re is also evidence from previous studies o f early-type 
galaxies (Helsdon et al. 2001; O'Sullivan et al. 2001a) that galax- 
ies in the centres of X-ray bright groups are affected by their en- 
vironment. They have a significantly steeper Lx'Lb relation, and 
are on average considerably more luminous than normal ellipticals. 
A large fraction of group dominant galaxies in one sample have 
been shown to have temperature profiles indicative of central cool- 
ing ( Helsdon et al. 2001), leading to the suggestion that their halos 
are actually the product of cooling flows associated with the sur- 
rounding group. Considering the differences between these domi- 
nant galaxies and their more normal counterparts, and the biasing 
effect their inclusion in samples of early-type galaxies seems to 
have, further investigation of the processes which have shaped their 
halos seems warranted. 

We have compiled a sample of 39 large. X-ray luminous 
early-type galaxies for which there is good quality ROSAT PSPC 
data available. We have analysed these data, and fitted two dimen- 
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sional, two component surface brightness profiles to them. We have 
also fitted two component spectral models, temperature, abundance 
and hardness profiles and produced three dimensional models of 
the galaxies. These allow us to model out contamination from sur- 
rounding cluster or group emission and the discrete source popu- 
lation within the galaxy. We can therefore examine the properties 
of the halo in detail, for the first time in a sample of this size. We 
can also compare the behaviour of this sample to that of samples 
of groups and clusters through relations between parameters such 
as temperature, optical and X-ray luminosity, velocity dispersion, 
surface brightness slope and gas entropy. In most cases, this is the 
first time these relations have been studied for halos at this mass 
scale. 

The paper is organised as follows. In Section |2| we describe 
our sample and the selection criteria used to create it. Section |3| 
gives details of the techniques used in reduction of the ROSAT 
PSPC data, and Section |4] describes the spectral and spatial fit- 
ting processes. Our results are presented in Section |5| with data 
from clusters and groups of galaxies included for comparison. We 
discuss the results and their implications in Section |S| and give 
our conclusions in Section Throughout the paper we assume 
Ho=50 km s~^ Mpc~^ , in order to simplify comparison with pre- 
vious studies of groups and clusters. Optical luminosities are nor- 
malised using the solar luminosity in the B band, Lbq = 5.2x 10'^^ 
ergs"\ 



2 SAMPLE SELECTION 

Our sample was selected from the Lyon-Meudon Extragalactic 
Data Archive (LEDA), specifically the PGC-ROM 1996 (2"^^ edi- 
tion). This contains information on ~ 100,000 galaxies, of which 
~40,000 have the necessary redshift and morphological data. 
Galaxies were selected to match the following selection criteria: 

(1) Absolute magnitude Ms < -19 

(2) Morphological T-type < -2 

(3) Virgocentric flow corrected recession velocity Vrec < 
10,500 km s"^ 

These criteria were chosen to produce a selection of optically lumi- 
nous nearby early-type galaxies. 

A list of galaxies matching these criteria was then compared to 
a catalogue of ROSATPSPC pointings, to produce an initial sample 
of galaxies with X-ray data. Only galaxies lying within the PSPC 
support structure (i.e. within ~30' of the pointing) were accepted, 
so as to ensure that the X-ray data were not strongly affected by 
vignetting effects or off-axis resolution problems. Pointings of less 
than 10 ksec were also ignored, as these were unlikely to provide 
X-ray data of sufficient quality. This initial sample contained 47 
galaxies. 

We then examined images of the raw X-ray data for each 
galaxy, to look for potential problems. In some cases we found that 
the galaxies appeared to be extremely compact or point-like, sug- 
gesting that surface brightness fitting would be difficult or impossi- 
ble. These objects were removed from the sample, as were galaxies 
in which an AGN or nearby quasar dominated the X-ray emission, 
to produce a final sample of 39 X-ray luminous early-type galaxies. 
It is worth noting that as the fraction of galaxies in which the halo 
was too compact or faint for analysis was small (~4%), it appears 
that the majority of massive early-type galaxies do possess bright, 
extended X-ray halos. TableQlists our targets. 



3 DATA REDUCTION 

Data reduction and analysis of the X-ray datasets were carried out 
using the ASTERIX software package. Before the datasets could be 
used, various sources of contamination had to be removed. Possible 
sources include charged particles and solar X-rays scattered into 
the telescope from the Earth's atmosphere. Onboard instrumenta- 
tion provides information which allows periods of high background 
to be identified. The master veto counter records the charged parti- 
cle flux, and we excluded all time periods during which the master 
veto rate exceeded 170 count s~^. Solar contamination causes a 
significant overall increase in the X-ray event rate. To remove this 
contamination we excluded all times during which the event rate 
deviated from the mean by more than 2a. This generally removed 
no more than a few percent of each dataset. 

After this cleaning process each dataset was binned into a 3- 
dimensional (x, y, energy) data cube. Spectra or images can be ex- 
tracted from such a cube by collapsing it along the axes. A model 
of the background was then generated based on an annulus taken 
from this cube. We used annuli of width 0. 1 ° , and inner radius 0.4° 
where possible. In cases where this would place the annulus close 
to the source we moved the annulus, generally to r = 0.55°. To en- 
sure that the background model was not biased by sources within 
the annulus, an iterative process was used to remove point sources 
of > 4.5 a significance. A number of our galaxies are found within 
groups and clusters of galaxies, many of which have their own X- 
ray halos. Our intention was to model these spectrally and spatially 
in order to accurately remove the effects of their contamination 
of our target galaxies. We therefore moved the annulus outward 
to avoid the emission, where possible. In cases where the emis- 
sion appeared to extend to the edge of the field of view, we used a 
background annulus at r=0.9°. This occurred for a small number of 
galaxies which lie in the centres of clusters (e.g. NGC 1399). The 
use of a background annulus which lies within the cluster emission 
means that we are likely to overestimate the true background and 
hence over-correct for it. However, as we are using the largest an- 
nulus possible, we should mimimize the degree of overestimation. 
We can also expect the central galaxy component of the emission 
to have a much higher surface brightness than the cluster emission, 
so that oversubtraction will have a negligible effect on it. Surface 
brightness fits should therefore be accurate for the central compo- 
nent, which is our main interest, and as good as is possible for the 
cluster component. 

The resulting background model was then used to produce a 
background- subtracted cube. Regions near the PSPC window sup- 
port structure were removed from these images, as objects in those 
areas would have been partially obscured during the observation. 
The cube was further corrected for dead time and vignetting ef- 
fects, and point sources were removed. 

Examination of background subtracted images allowed us to 
locate each galaxy and produce a radial profile of the surrounding 
region. From these profiles, regions of interest were selected, from 
which images or spectra for use in fitting could then be extracted. 
The majority of our galaxies are known to be members of groups 
or clusters, and as such we expected to see emission from an in- 
tergalactic medium (IGM) surrounding them. In cases where the 
galaxy did not appear to be contaminated with other emission, we 
defined the region of interest (Rol) as being within the radius at 
which the emission dropped to the background level. This region 
was suitable for both spatial and spectral fitting. In cases where 
contaminating intergalactic emission was seen, we defined sepa- 
rate regions of interest, one for spectral and one for spatial fits. For 
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Name 


RA 

(2000) 




DEC 
(2000) 


a 

(km s 


-1) 


(kms-i) 


D 
(Mpc) 


Re 

(') 


T 


Environment 


ESO 443-24 


13 


01 


01, 


.6 


-32 26 20 


279, 


.9 


4970.1 


99.4 


0, 


.388 


-3.1 


BGG 


IC 1459 


22 


57 


09, 


.5 


-36 27 37 


321, 


.4 


1522.0 


28.3 


0, 


.644 


-4.7 


BGG 


IC 4296 


13 


36 


38, 


.8 


-33 57 59 


341, 


.2 


3587.9 


71.8 


0, 


.953 


-4.8 


BGG 


IC 4765 


18 


47 


19, 


.0 


-63 19 49 


288, 


.4 


4344.9 


86.9 


0, 


.239 


-3.9 


BGG 


NGC 499 


01 


23 


11 


.5 


-1-33 27 36 


264, 


.2 


4482.7 


82.8 


0, 


.346 


-2.9 


BGG 


NGC 507 


01 


23 


40, 


.0 


-1-33 15 22 


295, 


.8 


5015.7 


100.3 


1, 


.285 


-3.3 


BGG 


NGC 533 


01 


25 


31 


.4 


-1-01 45 35 


250, 


.0 


5411.4 


95.5 


0, 


.792 


-4.7 


BGG 


NGC 720 


01 


53 


00, 


.4 


-13 44 21 


237, 


.7 


1622.0 


31.2 


0, 


.659 


-4.7 


BGG 


NGC 741 


01 


56 


20, 


.9 


-1-05 37 44 


288, 


.4 


5529.9 


91.6 


0, 


.869 


-4.8 


BGG 


NGC 1332 


03 


26 


17, 


.3 


-21 20 09 


328, 


.1 


1355.8 


29.5 


0, 


.467 


-2.9 


Group 


NGC 1380 


03 


36 


26, 


.9 


-34 58 33 


240, 


.4 


1617.6 


27.2 


0, 


.659 


-2.3 


Cluster 


NGC 1395 


03 


38 


29, 


.6 


-23 01 40 


241, 


.0 


1516.1 


30.8 


0, 


.757 


-4.8 


BGG 


NGC 1399 


03 


38 


28, 


.9 


-35 26 58 


329, 


.6 


1211.2 


27.2 


0, 


.706 


-4.5 


BCG 


NGC 1404 


03 


38 


51, 


.7 


-35 35 36 


212 


.3 


1701.9 


27.2 


0, 


.446 


-4.7 


Cluster 


NGC 1407 


03 


40 


12, 


.3 


-18 34 52 


279, 


.3 


1612.0 


30.9 


1, 


.199 


-4.5 


BGG 


NGC 1549 


04 


15 


45, 


.0 


-55 35 31 


203, 


.2 


932.5 


21.7 


0, 


.792 


-4.3 


Group 


NGC 1553 


04 


16 


10 


.3 


-55 46 51 


167, 


.5 


805.8 


21.7 


1, 


.094 


-2.3 


BGG 


NGC 2300 


07 


32 


19, 


.6 


-1-85 42 32 


263, 


.0 


2249.7 


41.5 


0, 


.524 


-3.4 


BGGt 


NGC 2832 


09 


19 


46, 


.5 


-1-33 45 02 


341, 


.2 


6992.2 


128.9 


0, 


.426 


-4.3 


BGG 


NGC 3091 


10 


00 


13, 


.8 


-19 38 14 


303, 


.4 


3670.2 


76.2 


0, 


.512 


-4.7 


BGG 


NGC 3607 


11 


16 


54, 


.1 


-1-18 03 12 


216, 


.8 


999.6 


29.7 


1, 


.094 


-3.1 


BGG 


NGC 3923 


11 


51 


02, 


.1 


-28 48 23 


269, 


.8 


1468.0 


26.8 


0, 


.889 


-4.6 


BGG 


NGC 4073 


12 


04 


26, 


.5 


-1-01 53 48 


267, 


.9 


5970.6 


119.1 


0, 


.931 


-4.1 


BGG 


NGC 4125 


12 


08 


07, 


.1 


-1-65 10 22 


239, 


.9 


1618.6 


38.9 


0, 


.998 


-4.8 


BGG 


NGC 4261 


12 


19 


22, 


.7 


-1-05 49 36 


316, 


.2 


2244.0 


47.2 


0, 


.644 


-4.8 


Cluster 


NGC 4291 


12 


20 


18, 


.1 


-1-75 22 21 


287, 


.7 


2043.9 


36.8 


0, 


.245 


-4.8 


Group 


NGC 4365 


12 


24 


27, 


.9 


-1-07 19 06 


268, 


.5 


1290.4 


23.9 


0, 


.830 


-4.8 


Cluster 


NGC 4472 


12 


29 


46, 


.5 


-1-07 59 58 


304, 


.8 


931.8 


23.9 


1 


.734 


-4.7 


BCG 


NGC 4552 


12 


35 


39, 


.9 


-1-12 33 25 


264, 


.2 


372.3 


23.9 


0, 


.500 


-4.6 


Cluster/AGN 


NGC 4636 


12 


42 


49, 


.8 


-1-02 41 17 


211, 


.3 


1125.2 


23.9 


1, 


.694 


-4.8 


Cluster/BGG 


NGC 4649 


12 


43 


40, 


.2 


-1-11 32 58 


342, 


.8 


1221.9 


23.9 


1, 


.227 


-4.6 


Cluster 


NGC 4697 


12 


48 


35, 


.9 


-05 48 02 


173, 


.4 


1232.2 


22.7 


1, 


.256 


-4.7 


BGG 


NGC 5128 


13 


25 


29, 


.0 


-43 01 00 


142, 


.6 


385.6 


5.8 


0, 


.708 


-2.1 


BGG/AGN 


NGC 5322 


13 


49 


15, 


.5 


-1-60 1 1 29 


239, 


.9 


2035.5 


41.7 


0, 


.587 


-4.8 


BGG 


NGC 5419 


14 


03 


38, 


.6 


-33 58 41 


329. 


.6 


4027.5 


80.5 


0, 


.723 


-4.2 


BGGt 


NGC 5846 


15 


06 


29, 


.3 


+01 36 25 


250. 


.0 


1890.0 


34.4 


1, 


.377 


-4.7 


BGG 


NGC 6269 


16 


57 


58, 


.4 


-1-27 51 19 


224, 


.4 


10435.0 


208.7 


0, 


.574 


-4.8 


BGG 


NGC 6482 


17 


51 


49, 


.0 


-1-23 04 20 


302, 


.0 


4102.0 


82.0 


0, 


.132 


-4.8 


Field ? 


NGC 7619 


23 


20 


14, 


.7 


-1-08 12 23 


310, 


.5 


3825.5 


60.0 


0, 


.536 


-4.7 


BCG 



Table 1. A list of galaxies included in our sample. RA and DEC are taken from the LEDA catalogue, as are morphological ty pe (T), Recession velocity (which 
is corrected for Virgocentric flow and movement within the local group), and velocity dispersion, a. Distances are taken from lPrugniel & Simie I] ll996l) where 
possible, or calculated from the recession velocity. Ho=50 is assumed in both ca ses, and for most objects, the d istance used is that of the group or cluster in 
which is resides. Optical effective radii. Re, are taken from Faber et al. ( 1989) and' Prugniel & HeraudeaJ il998^. Galaxy environment is taken from the group 
catalogues of Garcia 1 1993), Zabludoff & Mulchaev 1 1998) and White et al. 1 1999). BGG and BCG stand for Brightest Group Galaxy and Brightest Cluster 
Galaxy, indicating that the galaxy is the brightest (and presumably dominant) object in the surrounding system, t indicates galaxies which are not identified as 
brightest group galaxies in the group catalogues, but lie at the centre of the X-ray halo of their group. 



spatial fitting, we again define the Rol as being within the radius 
at which emission drops to the background level. This will contain 
both the galaxy and the surrounding group or cluster, allowing us 
to fit models to both and thereby accurately remove contaminat- 
ing emission. A background annulus for use with this region was 
selected as described above. 

For spectral fitting in cases where the galaxy is surrounded 
by contaminating group/cluster emission, a smaller Rol was de- 
fined, using the radius at which the galaxy emission dropped to the 
level of the surrounding group or cluster halo. Within this radius, 
emission should be dominated by components associated with the 
galaxy, though it may still be contaminated by group/cluster emis- 
sion along the line of sight. The emission outside this radius should 
be primarily produced by the group or cluster halo. We therefore 
take a local background spectrum from an annulus with an inner 



radius 0.05° larger than the new region of interest, and use this to 
generate a local background model. This local background should 
account for both the cosmic X~ray background in the region of 
interest, and for the group/cluster contamination along the line of 
sight. Depending on the form of the group/cluster halo, we might 
expect to under- subtract this contamination to some degree. For ex- 
ample, if the cluster halo is steeply declining in surface brightness 
outside our region of interest, the local background annulus will 
contain fewer counts and we would expect to underestimate the 
contamination along the line of sight. Ideally we would hope that 
any extended group or cluster halo would have a core radius some- 
what larger than the region of interest, so that its surface brightness 
is relatively constant over the whole area we are considering. Any 
serious under-subtraction of group or cluster component will have 
an effect on the spectral fits we obtain, particularly in the more 
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massive clusters where contamination by the surrounding hot ICM 
would produce fits with higher than expected temperatures. Simi- 
larly, if we have misjudged the radius at which emission associated 
with the galaxy becomes less important than that associated with 
the surrounding structure, we might expect to subtract part of the 
galaxy emission. Again, we would expect to see evidence of this in 
the results of spectral fits to the data. We return to this question in 
Section|5| 



4 SPECTRAL AND SPATIAL ANALYSIS 

Spectra for each galaxy were obtained by removing all data out- 
side the region of interest and collapsing the data cube along its 
X and y axes. As with the background annulus, an iterative pro- 
cess was used to remove point sources of >4.5 a significance in 
the region of interest, although any point sources within the D25 
diameter were assumed to be associated with the galaxy itself and 
therefore not removed. The spectra could then be fitted with a vari- 
ety of models. To provide a baseline for later fits and to measure the 
basic properties of the galax y halo, each spectrum was fitted with a 
MEK AL hot plasma model f Kaastra & Mewdll993| : lLiedahl et alj 
Il995h . Initially, only normalisation was fitted. Hydrogen absorp- 
tion colu mn densities wer e fixed at values determined from radio 
surveys liStark et alJll'993> . and temperature and metal abundance 
were fixed at 1 keV and 1 solar respectively. Parameters were then 
freed in order (temperature, hydrogen column, metallicity), and 
only re-frozen at their starting values if they became poorly defined 
or tended to extreme values. The basic temperature and metallicity 
values are likely to be representative of the majori ty of early-type 
galaxies (Matsushit a et alj200(ilMatsushitall200lh . but clearly fit- 
ted values are preferable. 

We then attempted to fit two component spectral mod- 
els for each galaxy. These generally included a power-law -1- 
MEKAL model and bremsstrahlung -1- MEKAL model in which 
the bremsstrahlung temperature was fixed at 7 keV. The first com- 
ponent was intended to represent a hard component produced 
by the population of X-ray binaries and other unresolved stellar 
sources within each galaxy. Recent Chandra studies of extragalac- 
tic X-ray binary populations suggest that this emission is reason- 
ably modeled by a power-law of index ~1.2 iSarazin et al. 200ll: 
iBlanton e t al. 2001), while ASCA studies show good fits using a 
high temperature bremsstrahlung model (Matsushita et al. 2000) . 
All models were fitted using the Cash statistic <CaslJll979l) .^e 
Cash statistic is defined as -21nL where L is the likelihood func- 
tion. This means that the most likely model has a minimum Cash 
statistic and that differences in the statistic are chi-squared (x^) dis- 
tributed. Thus confidence intervals can be calculated in the same 
was as for a conventional x'^ fit- By comparing the best fit Cash 
statistic for each model, and visually examining the spectral fit, we 
selected the best fit model for each galaxy. From this we could ex- 
tract (in most cases) the X-ray temperature and metallicity of the 
galaxy halo, as well as the X-ray flux from the galaxy halo and the 
stellar contribution. 

For each galaxy in the sample, we also derived simple pro- 
jected temperature and hardness profiles. Temperature profiles 
were produced by splitting the larger, surface brightness region 
of interest into several annuli, from which spectra were extracted. 
These spectra were then fitted using the best fitting spectral model 
for the galaxy as a whole. Initially the models were fitted with the 
metallicity and hydrogen column density frozen at their global best 
fit values. However, if the data quality permitted, we freed these 



parameters, providing us with crude metallicity profiles for a frac- 
tion of our sample. Given the limited spectral range of ROSAT, 
the inability of the PSPC to resolve individual spectral lines, and 
the small number of counts in each annulus, the abundances fitted 
should not be taken as accurate measurements. However, in some 
cases they do show interesting trends when considered in conjunc- 
tion with the temperature profiles. 

Hardness profiles were calculated in a somewhat similar man- 
ner. Again the larger region of interest was split in to a number of 
annuli. From each of these, counts in soft (0.3-1.3 keV) and hard 
(1.3-2.4 keV) bands were extracted and divided to produce a ra- 
tio of hard/soft emission. Simulated spectra indicate that a 0.5 keV 
MEKAL spectrum produces a value of ~ 0.5, while a power law 
of r=1.7 and a 7 keV bremsstrahlung spectrum produce values of 
~1.1 and 1.2 respectively. These profiles can be used to give a ba- 
sic idea of changes in emission across the galaxy and in particular 
to identify AGN. 

In order to study the spatial properties of the galaxy X-ray 
emission, we also performed fits to the 2-dimensional surface 
brightness profile of each galaxy. Following the initial data reduc- 
tion described in Section |3| we extracted an image in the 0.5-2 
keV band and corrected it for vigne tting. This was done using an 
energy-dependent exposure map (see lSnowden et alJl994l for a full 
description). Point sources were removed as in the spectral anal- 
ysis, and unrelated extended sources identified and excluded by 
hand. Use of the energy dependent exposure map results in a con- 
stant background level across the image, so a flat background was 
also determined and subtracted from the data. 

As in the case of spectral analysis, we can choose to fit a vari- 
ety of surface brightness models to our data. The most commonly 
used in this work was a modified King function (or "/3-profile") of 
the form: 

S{r) = Soil + (r/reo.e)')-"'^"+°-' (D 

where S{r) is the surface brightness at a given radius. So is 
the central surface brightness rcore is the core radius and /Sfu is a 
measure of the slope of the surface brightness profile. At various 
stages of the analysis we also fitted point source models and de 
Vaucouleurs r^^^ law models, using the form: 

S{r) = Se ■ exp{-7.67[(r/re)"''' - 1]} (2) 

where S^ is the surface brightness at re, the effective radius 
(the isophotal radius containing half the total luminosity.) 

Models were convolved with the PSPC point spread function 
at an energy determined from the mean photon energy of the emis- 
sion in the region of interest and then fitted to the data. Both spheri- 
cal and elliptical fits were possible when using the King and de Vau- 
couleurs models, with the position angle and major to minor axis 
ratio measuring the shape and orientation of elliptical fits. When 
using King models, all parameters (core radius, /Sfu normalisa- 
tion, X and y position and the ellipticity parameters) were usually 
allowed to vary freely, as were the parameters in point source mod- 
els (x and y position, normalisation). The de Vaucouleurs model is 
intended to represent the unresolved discrete source population of 
the galaxies, which we assume will take the same form as the stel- 
lar population. An alternative approach would be to assume that 
the discrete source population follows the distribution of globular 
clusters, but we lack accurate spatial models of this distribution for 
many of our galaxies. We therefore initially set the effective radius 
parameter of the de Vaucouleurs models to the value of the optical 
effective radius, and held it frozen in most cases. We did allow the 
effective radius to vary for a small number of galaxies, where the 
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Figure 2. Fitted metallicity profiles of NGC 1399 (left panel) and NGC 4636 (right panel) 



fit was well constrained by the data, in order to investigate differ- 
ences between the optical and X~ray stellar profiles. However only 
one galaxy, NGC 4697, was best fit by a model including a de Vau- 
couleurs component, demonstrating that in most of the galaxies in 
the sample, either emission from hot gas dominates or the X-ray bi- 
nary population does not follow the stellar population. A Chandra 
observation of NGC 4697 has shown it to have a relatively small 
gas halo, with much of its emission contributed by point sources 
i lSarazinetail20 01). so the success of the de Vaucouleurs compo- 
nent in this case could is perhaps unsurprising. The Chandra data 
for this galaxy are best fit by a surface brightness model which 
includes a de Vaucouleurs component and a King model whose pa- 
rameters are such that it is flat, providing a fairly constant contri- 
bution over the area studied. 

The use of 2-dimensional datasets to fit the surface bright- 
ness distribution can result in a low number of counts in many of 
the data bins. Under t hese conditions fitting performs poorly 
iNousek & Shu e"l989') so, as in the spectral analysis, maximum 
likelihood fitting based on the Cash statistic was used. However, 
the Cash statistic gives no indication of the absolute quality of the 
fit, only the quality relative to other fits. In order to gain some es- 
timate of the true fit quality, we used a Monte Carlo approach, in 
which the best fit 1- and 2-component model was used to gener- 
ate 1000 images of the groups, to which Poisson noise was added. 
These were then compared to the original image, the Cash statistic 



determined, and a Gaussian fitted to the resulting spread of val- 
ues. By comparing the actual Cash statistic to this distribution of 
values, we were able to determine the probability that the model 
could have produced the data. We were therefore able to identify 
cases where the 2-component fit was no more likely to reproduce 
the data than the 1 -component fit, and discard the 2-component fits 
for these galaxies. 



5 RESULTS 

5.1 Spectral and spatial fits 

Table|2| shows the results of our spectral fits. As mentioned previ- 
ously, metal abundances from ROSAT PSPC spectra are inherently 
unreliable, due to the relatively poor spectral resolution of the in- 
strument. This is reflected by the large errors on some of our fitted 
values, and by the fact that in some cases we had to hold metallic- 
ity frozen in order to secure a stable fit. The temperature values are 
more reliable, and give a mean temperature of 0.67±0.29 keV. 

As discussed in Section|3| a poor choice of local background 
for our targets could result in spectral fits biased by inclusion of 
group or cluster emission, or accidental subtraction of some of 
the galaxy emission. The clearest sign of this bias would be un- 
usually high or low fitted temperature, significantly different from 
those found in other studies. Four of our galaxies have temperatures 
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above 1 keV, and one (NGC 4697) has a temperature lower than 
0.3 keV. These are outside the range commonly considered typical 
for elliptical galaxies, so we compare the results for these galaxies 
to those in the literature. NGC 507 has a temperature marginally 
above 1 keV. Previo us ROSAT and ASCA studies have found sim- 
ilar temperatures I Kim & Fabbiano 1995^, iMatsumoto et a l.' 1997**) 
and metal abundances ( Buote 2000), an d more recent Chandra data 
also s upports a temperature of ~ 1 keV iForman et all200ll) . lBuotel 
i2002 ) fits a two temperature model to XMM-Newton EPIC data 
for NGC 1399, and recovers temperatures of ~1.5 and ~0.9 keV 
within 1', with both components approaching a temperature of 1.3- 
1.5 keV at 3-10'. Our value of ~1.2 keV is quite comparable to the 
cooler component, considering the region from which our spectrum 
was extracted. NGC 4073 has the most extreme temperature in Ta- 
ble|2l kT=1.6 keV. Analysis of XM M-Newton EPIC data for NGC 
4073 and its surrounding group I O' Sullivan et al l l2003l in prep.) 
suggest a temperature gradient within the stellar body of the galaxy, 
with projected temperatures rising from 1.4 to ~2 keV. We also find 
a high temperature for NGC 6269, kT=1.4±0.2 keV. No XMM- 
Newton or Chandra results are available in the literature for NGC 
6269, but previous analys es of the ROSAT data hav e found temper- 
atures of 1.3±0.15 keV iDahlem & Thierinj200(]l) and 1.36±0.07 
keV <MulchaevetalJll996l) . These are identical within the errors 
with our result. Lastly, NGC 4697 was one of the first elliptical 
galaxies to be observed with Chandra, which showed it t o possess 
a rela tively cool gas component with kT~0.29 keV I Sa razin et alJ 
1200 ih . This is a fairly good match to our measured temperature of 
kT=0.24±0.2 keV, and it should be noted that at least part of the 
difference between these results may arise from Chandra calibra- 
tion issues. In general, these comparisons suggest that our method 
of background selection and data analysis produces fairly accurate 
fits to the data. 

When calculating the luminosities of our targets, we are faced 
with the difficulty that while we have surface brightness models 
which should provide an accurate estimate of the total number 
of counts from the galaxy halo, we only have spectral informa- 
tion for a smaller central region. Ideally we would be able to si- 
multaneously fit spectral and spatial data, giving a true luminosity 
for each component (e. e. iLlovd-Davies et all2000l) . In practice we 
have chosen to scale up the MEKAL component of the spectral fit 
to the number of counts found for the surface brightness model. 
This allows us to calculate a gas luminosity for the galaxy compo- 
nent, but ignores the contribution from discrete sources. This lumi- 
nosity will therefore be an overestimate of the true gas luminosity 
associated with the galaxy. However, we expect large optically lu- 
minous galaxies such as our targets to be almost entirely dominated 
by gas emission. With a small number of exceptions, the spectral 
fits confirm this, suggesting that in most cases the overestimation 
is small. It is also notable that because the bremsstrahlung compo- 
nent peaks at a higher energy than the MEKAL, a given luminosity 
corresponds to a smaller number of bremsstrahlung counts (in the 
ROSAT band)than it would for a MEKAL model. This means that 
our overestimate of luminosity is reduced, as the number of counts 
associated with the bremsstrahlung component, and assumed in the 
scaling to be part of the MEKAL component, will produce only a 
small increase in gas luminosity. 

Temperature profiles for our sample of galaxies are shown in 
figureQ The extent of the profiles varies, owing to the relative qual- 
ity of data, length of exposure and g alaxy distance. Comparison 
of these profiles with those shown in Helsdon & Ponman 1 2000) 
shows them to be similar in most cases where the samples over- 
lap. We have classified the galaxy profiles into four groups; cooling 



cores (e.g. NGC 1399), hot cores (e.g. NGC 720), isothermal (e.g. 
NGC 4697) and those where the data quality prevents a judgment 
(e.g. NGC 4552, where the errors on the outermost bin are large 
enough to make it suspect). Cool and hot core galaxies are selected 
under the requirement that their central bin must be hotter or cooler 
than an average outer temperature by at least 20%, and that the er- 
rors in Tx must be smaller than this amount. It is also possible to 
see evidence of AGN activity in some of the profiles, particularly 
in the case of NGC 5128, where the central bin has a tempera- 
ture of ~5 keV and the rest of the galaxy < 1 keV. An interesting 
feature of some of the better defined profiles which show central 
cooling {e.g. NGC 507, NGC 1399, NGC 4636) is that the temper- 
ature rises with radius to a value above the apparent outer mean 
temperature, and the falls back to that mean, producing a tempera- 
ture peak at moderate radii. As all of these galaxies are embedded 
in larger group or cluster halos, this peak may mark the boundary 
of a group or cluster scale cooling flow, or the point of interaction 
between the galaxy and its environment. The observations avail- 
able for NGC 1399 and NGC 4636 contain very large numbers of 
counts, allowing us to include metallicity as a free parameter in the 
profile fitting. Metallicity profiles of these two galaxies are shown 
in Figure|2| Despite the large errors in some bins, it is notable that 
both metallicity profiles follow the same structure as seen in tem- 
perature; a central trough, rising to a peak at moderate radius, with 
an outer region of relatively low abundance. The peak tempera- 
ture and metallicity occur at approximately the same radius in both 
cases. In order to check that correlations between temperature and 
abundance in the fits were not biasing the results we modeled the 
fit space for the two bins on either side of the apparent break in the 
profile of NGC 4636, calculating fit statistics at a range of temper- 
atures and metallicities. Comparing the confidence regions for the 
two points shows that they are dissimilar to at least 9a significance, 
strongly suggesting the break in the profile is real. Confidence re- 
gions for the two bins are plotted in Figure|3| 

Table|3|lists the results of the surface brightness fitting for the 
galaxy halos of our target galaxies. In the majority of galaxies we 
obtain good quality fits, with relatively small errors on the core ra- 
dius and slope. As we are able to fit elliptical models, we also list 
the position angle of the major axis and axis ratio of the model fits. 
For five galaxies, we were unable to determine a reliable position 
angle, as the model was consistent (within errors) with being spher- 
ical. The best fit position angles of these galaxies are listed without 
errors, and should only be considered as rough estimates. 

5.2 The Lx-.Tx relation 

Figure|4]shows Lx plotted against temperature for our sample. As 
can be seen, we find a fairly tight relation between Lx and Tx, 
with onl jy a small num ber of outlier points. Using Kendall's K- 
statistic (Ponman '1982') to measure the strength of correlation, we 
find a significance o f ~4 .6a. We fitted the relation using ODR- 
PACK jBoggs et alJll989^ to perform orthogonal least squares re- 
gression, and found a slope of 4.8±0.7. This fitting method uses 
the errors in x and y for each point, but is unable to take of account 
of the fact that the errors are asymmetric. The mean of the upper 
and lower error is used for the error in each axis. 

This fitting method should be accurate as long as the points 
deviate from the mean relation only on account of the statistical 
errors. In cases where there is a large intrinsic scatter about the re- 
lation, the statistical errors do not provide an appropriate basis for 
the weighting of the data points. An alternate approach is to weight 
all the points equally, ignoring the statistical errors as misleading. 
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Table 2. Results of the spectral fits to our sample galaxies. Where possible, a absorbed MEKAL+bremsstrahlung (MK+BR) model was fitted, but in cases 
where the bremsstrahlung normalisation always tended to zero, this component was removed from the fit. All upper and lower ( 1 cr) errors on fitted parameters 
were calculated individually, but are shown as a single ± error when they are identical to two significant figures. Those galaxies for which metallicity could 
not be successfully fitted are listed with a fixed solar metallicity with no errors. Temperature profiles are classified as isothermal (1), cool core (C), hot core 
(H), or uncertain (?). 



To check our res ult -we also fit ted our Lx'Tx relation using the 
SLOPES package ilsobe et alll9 90) to perform an OLS bisector fit 
and ignoring statistical errors. The best fit slope for this technique 
was 5.05±0.44, very similar to the slope found when using the er- 
rors. 

From the temperature profiles for each galaxy, we are able 
to identify which of our targets have strong temperature gradi- 
ents which could affect the measured mean temperature. Excluding 
these 20 galaxies weakens the relation to -^2.80- significance, and 
gives a best fit slope (fitted using statistical errors) of 5.9±1.3. We 
also fit the sample of galaxies with known temperature gradients, 
and found a ~2.7cr correlation, with a slope of ^3.7. 

The Lx'-Tx relation has been used extensively in the study 
of groups and clusters of galaxies. Figure|5| shows our data points 



plotted alongside those for samples of groups faelsdon & Ponm^l 
mOQ) and cl usters ^avid et al. 1993; Mushotzkv & Scharl 199^ 
Ipairlev et a l. 2000). Our data follow a relation of similar slope to 
that of the groups (Helsdon & Ponman find a best fit slope of 
4.9±0.8), but offset to a lower luminosity or higher temperature. 
For a given temperature, our galaxies are a factor of ~3 less lumi- 
nous. Because of the scatter in both samples, there is some over- 
lap between the groups and galaxies, and some of our galaxy data 
points lie above the group best fit line. Conversely, the best fit re- 
lation for clusters is significantly shallower than that for groups or 
galaxies, though again there is a small region of overlap between 
the most luminous galaxies and the faintest clusters. 
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Name 



Core Radius 
(arcmin) 



Axis Ratio 



Position Angle 
(degrees) 



Rol Radius 
(degrees) 



Model 



ESO 443-24 


0.106 ±0.004 


0, 


IC 1459 


0.17 ±0.02 


0, 


IC 4296 


0.080 ± 0.004 


0, 


IC 4765 


0.39 ±0.02 


1, 


NGC 499 


0.85 ± 0.04 


0, 


NGC 507 


0.82 ± 0.03 


0, 


NGC 533 


0.058 ± 0.002 


0, 


NGC 720 


0.15 ± 0.02 


0, 


NGC 741 


0.172 ± 0.006 


0, 


NGC 1332 


0.010 ± 0.001 


0, 


NGC 1380 


0.090 ± 0.011 


0, 


NGC 1395 


0.23 ±0.03 


0, 


NGC 1399 


0.11 ±0.01 


0, 


NGC 1404 


0.32 ±0.04 


0, 


NGC 1407 


0.18 ±0.02 


0, 


NGC 1549 


0.021 ± 0.003 


0, 


NGC 1553 


0.43 ± 0.07 


0, 


NGC 2300 


0.23 ±0.02 


0, 


NGC 2832 


0.0100 ± 0.0003 


0, 


NGC 3091 


1.17 ±0.05 


1, 


NGC 3607 


0.63 ± 0.07 


0, 


NGC 3923 


0.010 ±0.001 


0, 


NGC 4073 


0.072 ± 0.002 


0, 


NGC 4125 


0.017 ±0.001 


0, 


NGC 4261 


0.37 ±0.03 


1, 


NGC 4291 


0.38 ±0.04 


0, 


NGC 4365 


0.54 ±0.08 


0, 


NGC 4472 


0.25 ±0.04 


0, 


NGC 4552 


0.098 ±0.014 


0, 


NGC 4636 


0.40 ± 0.06 


0, 


NGC 4649 


0.13 ±0.02 


0, 


NGC 4697 


0.99 ±0.15 


0, 


NGC 5128 


0.90 ±0.53 


0, 


NGC 5322 


0.0100 ± 0.0008 


0, 


NGC 5419 


4.5 ±0.2 


0, 


NGC 5846 


1.3 ±0.1 


0, 


NGC 6269 


1.21 ±0.02 


0, 


NGC 6482 


0.162 ±0.007 


0, 


NGC 7619 


0.031 ± 0.002 


0, 



0.55 ±0.03 
a+0.36 
"-0.09 
R+O.IS 
°-0.02 
d+0.6 

"'-o.oe 



0.54 ± 0.03 

q + 0.04 

■^-0.02 

on + 0.010 

"■^-0.009 
cr+0.64 

■^-o.oe 

-O+0.009 
^"-0.008 
1+0.04 
^-0.03 
9+0.06 
^-0.03 
9 ±0.02 
71-1+0.005 
'^-0.004 
R+0.02 
"-0.01 



.+0.06 
'-0.09 
1+0.12 
-0.06 

0.314 ±0.007 

^•°'J-0.04 



cc+O.OS 
^°-0.01 

.0+0.01 

^°-0.08 

,+1.4 

^-0.2 

tr7+0.04 

^"-0.03 

gQ + 0.04 

"'-'-0.03 
C-Q7+0.009 
-0.008 
gQ+0.02 
"^-0.03 
p.oc+0.007 
3 '''^-0.006 

567 ± 0.00 
46+0-13 

*"-0.07 



0.50 



+0.18 
0.07 
+0.05 
-0.04 
+0.04 
-0.02 
^+0.009 
^-0.001 
7+0.006 
' -0.005 



7+0.3 
-0.2 
H.4 
-0.2 
7+0.8 
-0.4 
1+0.1 
-0.7 

1 n5+°°3 

1.24 ±0.06 
1 34+0.09 
^•''*-0.08 
1.41+°-°® 



1.7^ 
1.6^ 
2.2j 

l.l^ 



1.1 

i.r 



-0.1 
B+0.3 
'-0.2 
,+0.2 
-0.3 

1.3 ±0.1 

1 r,o + 0.04 

^■^•^-0.03 

^■^°-0.01 

1.20 ±0.06 

1 56+°-^° 
^■''"-o.og 

1+0.2 
-0.1 
1+0.2 
-0.1 

1.4 ±0.1 
1.31 ±0.04 

-I 17+0.07 
^■^'-0.06 
a+0.4 
'-0.2 
-1+0.06 
-0.05 
1 629"'"'''^ 



1.4-' 

l.l^ 



1.. 

1.20^ 



0.003 
^■*-0.2 
1.3 ±0.1 
2.0 ±0.2 

1.08 ±0.02 

1 7+0.2 
^■'-0.1 

1 02+°°2 
^■'-'^-0.09 

1.18 ± 0.03 
1.8 ±0.3 
2.21 ±0.07 

1.44 ±0.06 
1.15 ± 0.03 
1 19+"-^ 

^■^='-0.08 
1 10+007 
^■^■^-0.05 
5+0.07 
-0.06 



-,+30 
-20 



306 ± 7 
276 ± 7 

190 
290j 
270 

240 ± 10 

294^7° 

12811° 

295 ±4 

210 

41+9 

176 ±3 

-10 
-6 
1+7 



-10 
-9 



1.18^ 



294^ 

91- „ 
174fl 

304^ 

330 

357 ±6 

325lg 

350 ± 10 

54l« 

265lg 

269lg 

37 ±5 

31 ±4 
83110 
164 ±5 
24 

26 ±6 
28l«, 
52.0 ±0.9 

44 ±6 

38 ±15 
220 ± 10 
308 ± 10 



0.045 
0.065 
0.12 
0.17 
0.25 
0.14 
0.16 
0.11 
0.12 
0.080 
0.045 
0.10 
0.12 
0.20 
0.20 
0.090 
0.17 
0.16 
0.15 
0.15 
0.14 
0.080 
0.10 
0.070 
0.10 
0.065 
0.090 
0.24 
0.045 
0.24 
0.080 
0.045 
0.25 
0.057 
0.20 
0.080 
0.10 
0.10 
0.20 



KI 
KI+KI 
KI+KI 
KI+KI 
KI+KI 
KI+KI 
KI+KI 
KI+KI 
KI+KI 
KI+KI 

KI 
KI+KI 
KI+KI 
KI+KI 
KI+KI 

KI 
KI+KI 
KI+KI 
KI+PS 
KI+KI 

KI 
KI+KI 
KI+KI 
KI+PS 
KI+KI 

KI 

KI 
KI+KI 

KI 
KI+KI 
KI+PS 
KI+DV 
KI+KI 

KI 
KI+PS 
KI+PS 
KI+PS 
KI+PS 
KI+KI 



Table 3. Results of the surface brightness fits to our galaxies, for the components associated with the galaxy halo. Best fit position angle values without errors 
are given for those galaxies where the angle was essentially unconstrained. All errors are quoted at the Ict. 



5.3 Lx'-Lb and Lb'-Tx relations 

One of the more common relations used in studies of early-type 
galsixies is the Lx'Lb relation. Numerous studi es based on ROSAT , 
ASCA or Einstein data have been published (e.^. Beuins et al. '199^; 
[Matsushita 2001; Brown & Bresman 1998; Fabbiano et al. 1992), 
and we have previously examined this relation in some detail in 
lO'Sullivan erS]j2001ah . to which we direct readers for a full dis- 
cussion of the relation and the effects of galaxy environment on 
it. Figure |S| shows the Lx'-Lb relation for our galaxies. For the 
sample as a whole, there is a 3.9a correlation, with a slope of 
~2.7. This is quite a steep relation, comparable to that foun d for 
a sample of BGGs in previous work JO'Sullivanetall2001iJ) . We 
have also pl otted the best fit relation fou nd for X-ray bright galaxy 
groups, from Helsd on & Potmia n ( 2002). The slope of this relation 
(2.6±0.4) is very similar to that found for our sample of galaxies. 
However, our relation is offset from that for groups, with galax- 



ies having X-ray luminosities a factor of ~9 higher than those of 
groups with equal optical luminosity, or conversely Lb values -^2.3 
times lower than groups of similar Lx. This result can be compared 
to the Lx'-Tx relation shown in Section l7!2l in which galaxies are 
offset to lower luminosities at a given temperature, compared to 
groups. 

The Lb'Tx relation for our galaxies is shown in Figure 
Once again, for this relation we find a fairly strong correlation 
(~4cr significance). The slope of the relation is comparable to 
that found for groups and clusters; 1.91±0.33 for our galaxy sam- 
ple, 1.64±0.23 for galaxy groups (Helsdon & Ponman 2002), and 
-^1.5 for galaxy clusters (iLlovd-Davies & Ponman .2002) . How- 
ever, where the relations for gro ups and clusters are essentially the 
same ( Helsdon & Ponman '2002), our relation for galaxies is signif- 
icantly offset to higher temperatures (by a factor of ~2, compared 
to groups) or lower Lb (by a factor of ~3). 
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Figure 4. Log Lx plotted against temperature for our sample of galaxies. Galaxies which show signs of central cooling are marked by open squares, those 
with central heating by stars, those which are apparently isothermal with circles, and those which are unclassified by triangles. The solid line shows the best 
fit line for the complete sample 



5.4 The a:Tx relation 

For undisturbed objects with hydrostatic halos, both X-ray temper- 
ature and velocity dispersion should be estimators of total system 
mass. It should be noted that the a quoted for our galaxies is a stel- 
lar velocity dispersion measured in the core of each object, whereas 
similar relations for groups and clusters use the velocity dispersions 
of the galaxies within those structures. Figure|8|shows the (j:Tx re- 
lation for our sample, again subdivided by temperature structure. 
For the sample as a whole, we find a ~3.1a correlation, with a 
slope of a ocTx"'^''*'' ''^. However, there appears to be a large de- 
gree of scatter about this line. Using the errors on the data points to 
measure the expected statistical scatter, we can estimate the intrin- 
sic scatter of the data. We find that the data points are ~1.8 times 
more scattered than would be expected from the statistical errors 
alone, hence the statistical scatter accounts for 56% of the variance 
we see. 



Also marked on the plot is the best fit relation for clusters of 
galaxies, taken from lWhite et alJ i 19991) . Our relation is indistin- 
guishable from this relation in both slop e and intercep t. The relation 
for galaxy groups is somewhat steeper: Eielsdon ifePon iuan 1 2000j) 



find a slope of cr oc Tx 



for their sample. Finally, Figure |8| 



shows a line representing /3spec=l, which is expected for systems 
where there is equipartition of specific energy between stars and 
gas. Our relation is consistent with this line, within the errors. 

Although there are many previous studies of samples of early- 
type galaxies, very few have measured temperatures for their tar- 
gets. Davis & White ( 1996), one of the few exceptions, fit a a:Tx 
relation to a sample of 26 galaxies observed using the ROSAT PSPC 
or Einstein IPC. They find a slightly steeper relation than ours, 
with a oc rx" '"'*" \ Adding the errors in quadrature we see that 
their measured slope is comparable with ours, at the la confidence 
level. Their best fit line is however offset to higher temperature. 
We believe this is probably caused by contamination by surround- 
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Figure 5. The Ly^.Ty^ relation for our sample, with similar relations for groups and clusters. Our galaxy data poin ts are marked (in black) by squares, and the 
relation by a solid line. Group and cluster data are shown in grey. Group data ai'e taken from Helsdon & Ponman 1 2000), and are marked by triangular points 
and a dashed line. Cluster data are drawn from David et al. ( 1993), Mushotzky & Scharl 1 1997) and Fairley et al. 1 2000), and are represented by circles and a 
dotted line. We find that the galaxy relation has a similar slope to that for groups, but is offset to fainter luminosities by a factor of ~3 



ing emission or discrete sources, as the spectral fits use a simple 
one component Raymond-Smith model (Raymond & Smith 1977). 
Unless the galaxies in their sample are completely dominated by 
halo emission, we would expect such contamination to raise the 
measured temperature of single component models, and we found 
evidence of such biases when comparing one and two component 
fits to galaxies in our own sample. 

5.5 /3/ it and Entropy 

Simple self-similar models of dark matter halos predict that X- 
ray emission from gas in the halo will always take the same 
form. Assuming the hot gas in the system is in hydrostatic equi- 
librium, the gas de nsity can be represented by a King profile 
iCavaliere & Fusco-Fe miano 1976). /3/it values for groups and 
clusters typically lie between 0.4 and 1 . Studies of low mass clus- 



ters and galaxy groups show that their surface brightness profiles 
become shallower with decreasing Tx iPonman et al. 1999). This 
is usually taken as an indication of additional physical processes 
affecting the gas, the result of which is that the the gas halo appears 
more extended and diffuse. Self-similar models also predict that 
gas entropy will vary linearly with system temperature, as entropy 
is here defined as S=Tx/ne''^ (where is the electron density 
of the plasma), and mean gas density will be constant for all sys- 
tems which virialised at the same redshift. For high mass clusters, 
this prediction matches the observed relation, but in lower mass 
systems the behaviour alters, with the trend flattening so that sys- 
tems of temperature ~1 keV seem to have entropy values (mea- 
sured at_^J_^tnriai]_scattered around a mean value of -^140 keV 
cm^ jLlovd-Davies et alj200ol) . 

The most common suggested processes whic h could affect 
the gas halo are heating, either by galaxy winds jPonman et alJ 
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Figure 6. Lx plotted against Lb for our sample of galaxies. Although there are some outlying points, a clear trend is visible. Galaxies in X-ray bright groups 
are marked by circles, those in X-ray faint groups by squares. The solid line marks our best fit relation, while the dotted line shows the best fit relation for 
X-ray bright groups iHelsdon & Ponman.200Z) . 



or AGN jw u et all2 00(i). or cooling of v ery low entropy gas 
/Muanwong et al."200l|).l Voit & Brvanl i200lh have recently sug- 
gested that a combination of cooling and star formation is likely to 
be a highly efficient way of producing the observed effects, as the 
heating will be centred in regions containing the lowest entropy gas, 
giving the maximum entropy increase. The existence of an entropy 
"floor", suggests that the level of entropy increase may be similar 
over a wide range of systems. In more massive systems, the en- 
tropy increase from shock heating is much larger than this amount, 
so it passes unnoticed. Only in small systems does it become the 
dominant contribution. Given this, and the fact that the suggested 
methods of raising the entropy are likely to occur predominantly 
within galaxies, we might expect that galaxy halos would show the 
effects of entropy increase very clearly. 

The entropy floor observed in low mass systems also has an 
effect on their surface brightness profiles. Raising the entropy of 
the intra-cluster medium (ICM) through heating will push gas out 



to larger radii if it occurs after the system forms, or prevent it from 
collapsing as far as expected if it occurs beforehand. The halo will 
therefore be more extended, with a lower central density, present- 
ing a flatter surface brightness profile. We therefore expect /3/,t to 
decrease with decreasing temperature, and this is observed across 
a wid e range of systems jPonman et alJll99^lLlovd-Davies et alJ 
|200(l). Given the likelihood that entropy increase has occurred in 
galaxy halos, we may also expect to find a relation between /3/it 
and temperature. 

Figure |9| shows the slope parameter Pfu plotted against tem- 
perature for our sample. There is no obvious trend in the points, 
and we find no statistically significant correlation. The scatter on 
the points is quite large, particularly at intermediate temperatures. 
There is no clear segregation of galaxies by temperature structure, 
and all classes show comparable amounts of scatter. There is some 
suggestion that AGN, BGGs and BCGs are more scattered than the 
more normal ellipticals, particularly if the normal galaxy with the 
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Figure 7. Lb plotted against temperature for our sample of galaxies. Circles mark galaxies which lie in X-ray bright gro ups or clusters, squares those in 
X-ray faint groups. The solid line is the best fit relation, while the dotted and dash-dotted lines show the relations for clusters iLlovd-Davies & Ponmari2002ft 
and X-ray bright groups tHelsdon & Ponman.20QZ) respectively. 



highest value of Pfu is excluded. This galaxy is NGC 1404, which 
is commonly considered to have suffered ram-pressure stripping of 
its halo. The sharp cut off in its surface brightness profile could con- 
ceivably have produced an unusually steep fit. However, the small 
size of the subsample makes any detailed comparison unreliable. 
The subsample seems to be centred on Pfu ~0.55, which is simi- 
lar to the mean value of our sample as a whole. 

Assuming that our galaxies are isothermal, we can extrap- 
olate from our 2-dimensional surface brightness models to 3- 
dimensional density models of our galaxies. Equation Q describes 
the 2D models and we can describe the 3D models as shown in 
Equation|3| 

. 3/3 



P(r) = P(0)[1 + (r/rcore) ]' 



(3) 



The density normalisation, p(o)' can then be determined from the 
surface brightness normalisation, assuming the temperature and 
metallicity determined from spectral fitting. From these 3D mod- 



els, it is possible to derive gas properties such as density, entropy, 
cooling time, etc, as a function of radius. In order to be able to 
compare the resulting profiles fairly, we need to view them on a 
common scale. This can be done by scaling the profiles by the 
virial radius of the system, or more usually to fixed overdensity 
radius, such as R2oo- The overdensity is calculated relative to the 
critical density of the universe at the redshift of formation. This is 
unknown for most systems, and in the case of clusters and groups 
is generally taken to be the redshift of observation. For galaxies, 
such an assumption is very unlikely to be accurate, even taking 
into account the fact that early-type galaxies will have potentials 
determined by the density of the universe at the redshift of last ma- 
jor merger rather than that at which the majority of their stellar 
population formed. In order to calculate R200 , w e assume a mean 
redshift of formation for ellipticals of 2/orm=2 iKauffmann et alj 
1996; van Dokkum & Franx 2001). Using this v alue, we can cal- 
culate R200 as described in IBalogh et alJ h999h and iBabul et"ai] 
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Figure 8. Velocity dispersion, cr, plotted against temperature for our sample of galaxies. Galaxies which show signs of central cooling are marked with open 
circles, those with central heating by stars, relatively isothermal galaxies by squares and those we have been unable to classify by triangles. The best fit line to 
the sample as a whole is marked by a so lid line, and the fit s to the two subsamples are marked by dashed lines. /3spec=l is marked by a dotted line, while the 
best fit relation for galaxy clusters (from lWhite et all 19991) is marked by a dot-dashed line. 



2002 ), taking variation of overdensity with redshift from lEke et alJ 

Figure ^| shows entropy, calculated at one tenth of R200 • 
Once again, galaxies are marked to show their different tempera- 
ture structures, but while the data shows a large amount of scatter, 
there is no significant correlation with temperature, or segregation 
by temperature structure. The galaxy data points do mainly fall at 
the low temperature end of the S :Tx trend for groups and clusters, 
but show no evidence of any trend themselves. 



6 DISCUSSION 

The results presented in the previous section do not lend themselves 
to a simple explanation. The halos of massive early-type galaxies 
seem to have some obvious similarities to those of galaxy groups 



and clusters, but also some intriguing differences. Clearly a num- 
ber of the relations which hold for groups and clusters also apply to 
early-type galaxies, the most important examples being the a:Tx 
and Lx'-Tx relations. The a:Tx relation is particularly striking, in 
that it agrees closely with the relation found for galaxy clusters. 
The Lx'Tx relation has a similar slope to that of galaxy groups, 
but is offset to lower luminosities. The L_g:Tx and Lx'-Lb rela- 
tions also show offsets from the cluster and group relations, but to 
higher X-ray luminosity and temperature. Another important result 
is the roughly constant value of f3fit across a wide range of tem- 
peratures. Although there is some suggestion of influence by sur- 
rounding groups and clusters on this parameter, the environment 
seems to increase scatter rather than producing a trend. Together 
with the results of our entropy calculations for the sample, the be- 
haviour of f3fit suggests that if preheating effects are important in 
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Figure 9. Plots of 13 against temperature for our sample of galaxies. On the left, galaxies are grouped by their temperature properties, with those with cool 
cores marked by open circles, those with hot cores by stars, those with isothermal profiles by squares, and those we have been unable to classify by triangles. 
On the right we segregate galaxies depending on their environment, with filled triangles representing brightest group galaxies, open triangles brightest cluster 
galaxies, stars AGN, and squares normal early-type galaxies. Error bars are omitted for clarity. 



1.5 I r 




0.5 1 1.5 

kT (keV) 

Figure 3. Confidence regions for bins 6 (solid lines) and 7 (dotted lines) of 
the temperature and metallicity profiles of NGC 4636. The contours show 
confidence levels of 1, 3, 6, 9 and I2a, and only at the I2a level do the two 
regions come close to overlap. 



galaxy halos, they produce results quite different to those seen in 
groups. 



6.1 Comparison with Groups and Clusters 

One of the most important results to arise from studies of galaxy 
groups and clusters is the demonstration that simple self-similar 
models do not describe low mass systems well. There is good evi- 
dence that galaxy groups differ from higher mass clusters in a num- 
ber of important ways. The most widely accepted explanation for 
these differences is that the hot gas in the halos of these systems 
is affected not only by the processes involved in formation of the 
system as a whole, but also by processes such as star formation, 
AGN heating, and gas cooling. A clear sign of this is seen in the 
behaviour of gas entropy in systems of different mass. In high mass 
clusters, entropy can be fairly accurately predicted from simple 
models in which the gas is heated (and has its entropy profile set) 
during formation of the system. As the cluster builds up, gas flows 
into the potential well and is shock heated to a degree dependent on 
the depth of the potential. This dependence of halo temperature on 
system mass can be seen in the criTx relation for clusters, in which 
Tx is strongly correlated with a, with Tx °5 (^^^ as expected from 
arguments based on the virial theorem. 

However, as system mass (and temperature) drop, the gas 
entropy begins to depart from the predicted relation, with low 
mass galaxy group s having higher than expected gas entropies. 
iLlovd-Davies et au llOOQ) have shown that at temperatures below 
~2 keV, gas entropy appears to remain roughly constant, scatter- 
ing about a mean value of 140 keV cm^. An effect is also seen at 
around this temperat ure in the surface brig htness profiles of galaxy 
groups and clusters jPonman et alJll99^ . High temperature sys- 
tems (Tx > 4 keV) have profiles which scale self-similarly, but as 
Tx is lowered, the surface brightness profiles observed are found to 
be shallower, with lower central densities. This has an obvious ef- 
fect on the measured X-ray luminosities of these systems, making 
them fainter than predicted. 
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Figure 10. Mean gas entropy measured at one tenth of the virial radius plotted against mean halo temperature. Galaxies with cool cores are marked as circles, 
galaxies with relatively isothermal halos by squares, those with hot cores by stars and those we were unable to classify by triangles. Grey points with eiTor 
bars represent groups and clusters taken from the Birmingham-CfA cluster scaling project (sec iSanderson et alj2Q02.) . The approximate level of the entropy 
floor for groups and clusters is marked by a dotted line. 



Both of these trends are explained through the effects of non- 
gravitational processes on the gas. To take a simple model, we 
can imagine a system at the time of formation. As stated above, 
we expect the gas falling into the potential of the system to be 
heated by gravitational processes, but we might also expect heating 
from other sources. For example, as galaxies form in the system 
we could expect star formation, galaxy winds, and AGN to affect 
the gas in the system halo. In a large system, the contribution from 
these sources will be small compared to that of gravitational heat- 
ing. However, in smaller systems it will be more important, and 
will eventually dominate the energy of the gas halo. This additional 
heating will have a number of effects, such as raising the gas tem- 
perature, causing the halo to expand by moving gas to higher radii, 
raising the entropy of the gas, and so on. In practice, it is possible 
to achieve these effects through a number of proces ses, or a com- 
bination of several. One promising model is that of IVoit & BrvanI 



in which low entropy gas cools rapidly, providing material 
for star formation. This star formation then not only removes low 
entropy gas from the system (raising the mean entropy), but pro- 
vides heating which is focused in the areas in which low entropy 
material dominates. From our point of view however, the method is 
not important; all the suggested processes would occur in galaxies, 
and might be expected to occur preferentially in large galaxies at 
the bottom of a group or cluster potential well. For simplicity, we 
will refer to this as the preheating model. 

6.1.1 The CT.'Tx relation 

One of the clearest similarities between the relations for our galax- 
ies and those of larger structures is the correspondence of the aiTx 
relation to that of galaxy clusters. As discussed above, the cor- 
relation in clusters is expected, and shows that both Tx and a 
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are good measures of the potential. In galaxy groups, the slope 
of the relation is observed to be somewhat steeper, ~1.7±0.3 
iHelsdon & Pon man 200(|). One explanation for this steeper slope 
is that Tx is raised in low mass systems by the heating and re- 
moval of cool gas described above. An alternative is the suggestion 
that the velocity dispersions measured in groups could be biased. 
There are a number of groups which, despite apparently possessing 
extended X-ray halos, hav e exceptionally low velocity dispersions 
iHelsdon & Ponmarl200(t) . If these a values were accurate, the po- 
tential of the group would not only be too shallow to produce the 
halo luminosity observed, but would be too shallow for the group 
to have collapsed within the age of the universe. Several reasons 
for underestimation of a can be postulated. It is possible that these 
low mass groups form a central core of bright galaxies, with fainter 
members at higher radii. The brightest galaxies are most likely to 
be recognised as group members and have measured redshifts, so 
they will dominate any calculation of a. Tidal interactions between 
group members could reduce velocity dispersion by transferring 
orbital energy from the galaxies to their stars. It has also been sug- 
gested that low mass groups may be biased because the groups of- 
ten have prolate structures. For groups whose major axis lies near 
the plane of the s ky, this particular s tructure could lead to an under- 
estimation of a iTovmassian et aljl20 021. biasing samples includ- 
ing such low temperature systems. If a is biased, correcting the 
measurements for the lowest temperature groups would probably 
shift the best fit relation into agreement with that seen in clusters. It 
is worth noting that the groups with the lowest velocity dispersions 
listed injilelsdon & Ponman (2000) are also amongst the poorest; 
all groups in their sample with a < 150 km s~^ have only 3-5 
member galaxies. Calculating velocity dispersion from such small 
numbers can introduc e a bias, produci ng values which are underes- 
timated by up to 15% jHelsdorl2o'o3) . 

The fact that our results agree well with the cluster relation 
suggests that again, we are looking at systems in which both Tx 
and a are good measures of the potential. From a preheating point 
of view this is surprising. If we expect preheating processes to oc- 
cur in large galaxies, then we might expect the Tx of their halos 
to be raised, in much the same way as we see in galaxy groups, 
producing a steeper relation. We also need an explanation of why 
the gas in these galaxies has a temperature which is related to the 
depth of their potential well. As well as additional heating from star 
formation and AGN activity, we expect large quantities of gas to be 
lost from the stars in the galaxy. As this gas is produced within the 
galaxy, we cannot expect shock heating during infall, so we might 
initially believe its temperature to be entirely determined by super- 
nova heating. 

iHelsdon et alj ^2001) consider a related question, that of the 
relative importance of different energy sources which contribute to 
the total X-ray luminosity of early-type galaxies. Their Figure 8 is 
aplot of Lx/Lb against Lb, showing the contributions to Lx from 
discrete sources, SNIa and gravity. From our point of view the con- 
tribution from discrete sources is irrelevant, as we are interested in 
energy input to the gas in the galaxy. The gravitational contribu- 
tion is a combination of two processes, firstly a contribution from 
the velocity of the stars in the potential (gas lost from these stars 
will have an added kinetic energy component from their velocity, 
which will be thermalized in the surrounding ISM), and secondly a 
contribution from work done on the gas as gravity causes it to con- 
tract and cool. The important result with regards to our situation is 
that while the SNIa energy input scales with Lb, the gravitational 
input scales with a'^. This means that for low mass systems, the 
dominant contribution to Lx is from supernova heating, but above 



Ls ~ 10 1/0 gravitational work begins to dominate. If we con- 
sider these two factors as energy inputs to the gas rather than as 
contributors to Lx, we can see that the gas temperature is likely to 
be determined by the supernova rate in low mass systems, and by 
the depth of the potential in high mass systems. The point at which 
the two are equal will depend on the details of the model, but if 
we follow the assumptions made by Helsdon et al. 1 2001), all our 
galaxies lie in the high mass, gravitationally dominated region of 
the plot. We could therefore expect gas temperature to depend 
on the depth of the potential, whether the gas has an internal 
or external origin. 

It seems likely, from this result, that we can draw similar con- 
clusions from the a:Tx relation in galaxies as we would in clusters. 
Tx and a are both probably good measures of mass. Preheating, by 
whatever method, does not appear to have the effect on galaxy halos 
that it has on those of groups. Like clusters, the galaxy relation is 
consistent with the systems having Pspec = 1, suggesting that there 
is equipartition of specific energy between stars and gas. As we will 
discuss later, this has important implications, in that it suggests that 
the optical and X-ray density profiles should be similar. One fur- 
ther consideration is the intrinsic scatter in the points. Our data has 
-^1.8 times as much scatter as we would expect from statistical er- 
rors alone, giving us a non-statistical unceratinty in T of ±0.04 for 
any given value of a. For comparison, the galaxy groups studied 
by Helsdon & Ponman (2000) are less scattered, having only ~1.4 
times as much as would be expected from the errors, sothe non- 
statistical uncertainty in T for these systems is ±0.86. The degree of 
scatter in the galaxy data could have a number of causes. Gas tem- 
perature could be affected by many processes, related to the galaxy 
or the surrounding environment. Velocity dispersion could also be 
affected by processes associated with the formation or merger his- 
tory of the galaxy. However, it would appear that galaxy groups 
have a larger scatter in properties, suggesting that the galaxies are 
less affected by external influences. 

6.1.2 The Lx.'Tx relation and (3fit 

Accepting a and Tx as indicators of the mass of the system, we 
next consider the Lx '.Tx relation. Here, we find that the slope of the 
relation is steeper than that of clusters, as steep as that of groups. 
The relation is also offset, so that at a given temperature, galaxies 
have a lower luminosity than groups. The steep slope is usually 
explained as a product of preheating - the additional heating of the 
gas raises the temperature slightly and moves gas to higher radii, 
lowering the central density and therefore Lx. If we assume that the 
steep slope seen in galaxies is caused by the same processes which 
cause it in groups, then this relation is a strong piece of evidence 
for the effects of preheating in galaxies. 

However, ther e are other ways in w hich we might produce 
such a steep slope. 'Helsd on et alj i200ll) found a strong correla- 
tion between the X-ray properties of group dominant galaxies and 
those of the groups in which they are found. Given the signs of cen- 
tral cooling in many of these groups, they suggested that what had 
been initially identified as the halos of the dominant galaxies were 
in fact group scale cooling flows, centred on the dominant galaxy 
because it lies at the bottom of the group potential. They also found 
that the Lx of this central galaxy halo/cooling flow was ~25% of 
the Lx of the group. As 29 of our 39 galaxies are dominant galax- 
ies in groups, clusters or cluster subclumps, we must consider the 
idea that our relations could be dominated by cooling flows. In that 
case, we would expect the Lx'.Tx relation to have a similar slope 
to that of groups, but to be offset to lower Lx values by a factor of 
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4, and also to lower temperatures. Provided the temperature drop 
is not too large, this could reproduce the relation we see very well. 
On the other hand, we do not see any segregation in the data be- 
tween those galaxies which show signs of central cooling and those 
which do not, nor do we see any difference between galaxies at 
the centres of groups and those in other environments. This argues 
against cooling flows as the driver of the relation. We will discuss 
the evidence for and against cooling flows as the dominant factor 
in Section l6^ 

The lack of a relation between /3/it and Tx argues against both 
cooling flows and preheating as the source of the Lx'-Tx relation. 
In galaxy groups, preheating causes gas to move out to high radii, 
reducing Pfu as the surface brightness profile becomes flatter. As 
preheating is more effective in smaller mass systems, groups and 
low mass clusters show a correlation between (3 fa and Tx with 
cooler systems having flatter profiles. In galaxies, despite a very 
large scatter, we see no trend with temperature. Our galaxies lie 
around a mean value of /Sfu = 0.55, which means that even if the 
group and cluster /3/it:?x relation levelled off at low temperature, 
our sample would not be consistent with it. This strongly suggests 
that preheating is not the cause of the steep slope of the Lx'-Tx 
relation. 

If a change in the slope of the surface brightness profile is not 
responsible for the drop in Lx needed to produce the Lx'-Tx rela- 
tion, then there must be a drop in normalisation. This is demon- 
strated in Figure 1111 which shows the surface brightness pro- 
files of our sample, scaled so that if they were behaving self- 
similarly, they would coincide. Details of this scaling are given 
in the figure caption. A comparison between Figure fTTI and Fig- 
ure 1 of lPonman et alj il999.^ demonstrates the difference between 
the behaviour of our galaxies and that o f groups and clusters. The 
high temperature clusters in Figure 1 of 'Pon man et alj ^99^ be- 
have self-similarly, having roughly equal surface brightness slopes 
and scaled normalisations. The lower temperature groups move 
away from self-similarity, as their slopes flatten and central sur- 
face brightness (and density) is lowered. In the galaxies, although 
there is a large amount of scatter, the slope remains roughly con- 
stant across the range of temperatures, but the normalisation of the 
profiles seems to drop with decreasing temperatures. This suggests 
that whereas in galaxy groups the steepening of the Lx'Tx re- 
lation is caused by energy injection and movement of gas to 
higher radii, in galaxies is caused by changes in normalisation 
and therefore in the overall gas fraction of the systems. The 
processes responsible for this change are not clear, but this result 
demonstrates how differently galaxies behave compared to groups 
and clusters. 

It is interesting to note that the mean slope of the surface 
brightness profiles (/^/^ = 0.55) is similar to the mean slope of 
the optical surface brightness profile. Elliptical galaxies are well 
described in the optical by a de Vaucouleurs profile, which is sim- 
ilar (outside the central regions) to a King model with Pfu = 0.5. 
As mentioned in Section l6. 1.11 the a:Tx relation is consistent with 
Papec = 1, which leads us to expect a similarity between the op- 
tical and X-ray density profiles. The similarity in surface bright- 
ness profiles is therefore further confirmation of the a:Tx result. 
A further interesting note is that galaxy wind models, in which the 
majority of the gas in the halo is produced by stellar mass loss, pre- 
dict X-ray surface brightness profiles similar to those in the optical 
jPell egrini & Ciotti 1998). The exact value of f3fit would depend 
on the wind state of the system, with only relatively hydrostatic ha- 
los having Pfu = 0.5. Galaxies dominated by supersonic outflows, 
or those in which inflows are important, would be expected to have 



steeper profiles owing to their high central gas densities. Galaxies 
dominated by subsonic outflows would be expected to have slightly 
flatter profiles (/3/it < 0.5), as the outflowing gas has a higher den- 
sity at large radii than in a supersonic flow, and less of a central 
peak. If such models are applicable, it seems likely from the high 
mass of our galaxies that they would be in the inflow stage, so we 
would expect /3/it > 0.5, in agreement with our observations. 



6.L3 Entropy 

The results of our entropy measurements are in broad agreement 
with those from the surface brightness profiles. In galaxy clusters 
and groups we see evidence of similarity breaking and preheating, 
leading to a trend in entropy which levels off at the entropy floor. 
The galaxy data points do not show any sign of a trend, and al- 
though they may appear to be consistent with the general trend 
in higher mass systems, the scatter in the data is very large. Our 
method of calculating entropy has two important sources of scat- 
ter associated with it. Firstly, we must assume that our galaxies 
are isothermal, despite the fact that we know that many of them 
have temperature gradients. Secondly, we measure the entropy at 
one tenth of the virial radius, and when calculating the virial radius 
we must assume a redshift of formation (or last major merger). Al- 
though the assumed value of Zform = 2 is probably a reasonable 
mean, there will clearly be variation between individual galaxies. 
Our virial radii will therefore be inaccurate to some degree, mean- 
ing that we are actually measuring entropy at a range of scaled radii. 

Despite these difficulties, we can draw some important results 
from the entropy values. A large proportion of the data points lie 
below 100 keV cm^, meaning that they are below the entropy floor 
observed in galaxy groups. One possible reason for these low val- 
ues is that the galaxies form earlier than the surrounding groups. 
The density of systems virialising at a given epoch is related to the 
critical density of the universe at that time, and so the earlier an 
object virialises, the denser we expect it to be. As entropy is in- 

2/3 

versely proportional to , an equal amount of energy injected 
into a denser system will produce a smaller increase in entropy. 
Galaxies, forming at z ~ 2, might therefore be expected to have 
higher densities and a lower entropy floor than groups and clusters 
which have formed more recently. However, we might also expect 
that as the site of the processes responsible for entropy increase, 
the amount of energy available to affect entropy in galaxies would 
be larger than in groups. It is also worth considering that if galax- 
ies had formed with halos of very low entropy gas, it would have 
cooled on timescales considerably shorter than the Hubble time, 
probably leading to star formation, heating and a rise in observed 
entropy ( Vo it & Brvaa^OOD . 

An alternative viewpoint is that the majority of the gas in 
these systems is being produced by stellar mass loss rather than 
infall, so we should not expect entropy to behave as it does in 
groups. In this case, galaxy wind models should give a reasonable 
approximation of the behaviour of the halo, in the absence of sig- 
nificant environmental influences. Several published simulations of 
galaxy halo development give gas density and temperature profiles 
(Ciotti et al. 1991; Pellegrini & Ciotti 1998; Brighenti & Mathew^ 
1999) for their model galaxies. In models of outflowing winds, both 
gas density and temperature fall with radius, but density falls more 
rapidly, so we would expect an entropy profile which increases with 
radius. In models which have developed gas inflow (cooling flows), 
temperature may rise with radius throughout most of the model, so 
again we expect entropy to rise with radius. This matches what we 
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Figure 11. Scaled X-ray surface brightness profiles for our galaxies, arranged to show departures from self-similarity. The profile s have been scaled t o the 
virial radii of the galaxies (assuming a redshift of formation Zfaj.^=2), and by Tx^/^ to account for self-similar scaling with mass IPonman et alll99 9^. We 
have also converted from counts to intensity, using the best fit MEKAL component of each galaxy to allow for the temperature and metallicity dependence 
of the plasma emissivity. As all our galaxies are at low redshift (2 < 0.035) we ignore cosmological effects as neghgible. In each plot bold, dark lines show 
the galaxies in the energy band in question, naiTow, grey lines show those in the band above, and dotted lines those in the band below. Although there is some 
scatter, it is clear that all have similar slopes and that cooler systems have lower normalisations. 



observe in our measured entropy profiles. The entropy predicted at 
any given radius depends on tlie details of tlie model, e.g. system 
mass, age, supernova rate, mass injection rate, etc. Entropies rang- 
ing from 20-300 keV cm^ might be expected for galaxies such as 
those in our sample. The majority of our galaxies do have entropies 
within these limits, and considering the expected scatter, the agree- 
ment between models and measurements is fairly good. 

6.2 Cooling Flows 

A number of our galaxies have temperature profiles indicative of 
central cooling. This is not surprising, as they are fairly massive ob- 
jects, with large halos, and reside within larger structures which are 
themselves probably capable of producing cooling flows. Among 



the relations we have examined, several show behaviour which 
could be explained easily as the product of group cooling flows. 
The best example of this is the Lx:Tx relation, where the slope 
is identical (within errors) with that found for groups, but offset to 
lower X-ray luminosities. As discussed in Section lS. 1.21 this offset 
(a factor of 3 in Lx) is very similar to that we would pr edict (a 
factor of 4), based on studies of X-ray bright groups I Helsdon et alj 
I2OOII) . We would expect the measured Tx of the cooling flow re- 
gion to be lower than the group temperature and this could explain 
the difference between predicted and measured offset. We might 
also be able to explain the offsets observed in the Lx'Lb and 
Lb '.Tx relation using this model. In both cases we are comparing 
a parameter which is determined by the galaxy and only weakly 
influenced by the group (Lb) to a parameter which is determined 
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by the cooling flow and hence by the properties of the group (Lx 
or Tx). In such a situation, we would have to expect Lx and Tx 
to appear unusually high for the central galaxy, as they would actu- 
ally be related to the much larger system of the surrounding group 
or cluster. Comparing with optical luminosity, we would only ex- 
pect to find a relation without an offset if we compared them to 
a value of Ls calculated for all the galaxies in the group, rather 
than just the central dominant elliptical. We would also expect the 
galaxies in X-ray faint groups to behave differently, as their halos 
could not be the product of group scale cooling flows. Unfortu- 
nately, the small numbers of such objects in our sample means we 
cannot test whether they follow a different trend to those in X-ray 
bright groups, but they do indeed appear to fall at lower Lx and 
Tx , f or a given L b ■ 

However, there are two powerful arguments against group 
scale cooling flows as the dominating influence in our sample. 
Firstly, we find that a significant subset (one third) of our sample 
do not show signs of central cooling. They are instead approxi- 
mately isothermal, or have a central rise in temperature. This does 
not necessarily mean that they have not been at the centre of a group 
cooling flow in the past, but such a flow would have to have been 
disrupted, and so would be unlikely to produce the emission we 
assume to be a halo. Secondly, and most importantly, we see no 
sign of a difference between those galaxies with or without signs 
of central cooling. The only segregation we see in any of the re- 
lations is a tendency for galaxies which seem to harbour cooling 
flows to be higher mass systems, and so to have higher values of 
(T, Tx, I/x, etc. We see no other difference between galaxies with 
differing temperature profiles, and we specifically do not see the 
galaxies with apparent cooling flows driving the Lx'-Tx relation. 
This strongly suggests that while a number of our galaxies do har- 
bour cooling flows, some of them quite large, their halos are not 
simply cooling flows formed by surrounding groups and clusters, 
unless the properties of the gas halo are independent of whether it 
arises from stellar winds or group inflow. 

6.3 Stellar Mass Loss and Galaxy Winds 

Large scale cooling flows, in which most of the gas in the galaxy 
halo has an external origin, do not appear to provide a good model 
for our galaxy sample. The alternative is a model in which the ma- 
jority of the gas is produced (and heated) internally, via stellar mass 
loss. We have already discussed in Section lS.l.ll the relative contri- 
butions to the energy of gas produced in a galaxy from supemovae 
and gravitational processes. Our sample is made up of galaxies with 
high stellar masses, leading us to believe that the gravitational po- 
tential should be the dominant source of energy. This provides an 
important link between the X-ray properties of the galaxies and 
their mass, which could explain the correlations we see between 
optical and X-ray properties in our sample. If gravity is the dom- 
inant source of energy in the halos of our galaxies, we would ex- 
pect them to behave like clusters on many of the relations we have 
examined. The best example of this is the a:Tx relation, where the 
galaxies and clusters have best fit lines which are indistinguishable. 

However, there are a number of ways in which the relations 
do not behave like those of galaxy clusters, and we need to find ex- 
planations for these differences. The most interesting, and perhaps 
the most important, is the Lx'-Tx relation. This behaves like that 
of galaxy groups, rather than galaxy clusters. However, instead of 
a decrease in central gas density in low mass systems, leading to 
a trend in f3fit with temperature, we see a constant value of /Sfu 
and a decrease in the normalisation of the surface brightness pro- 



files at low masses. This suggests that the lower mass members of 
our sample lie progressively further below the cluster Lx'-Tx rela- 
tion because they have a lower gas fraction than their higher mass 
counterparts. To explain the Lx'Tx relation we observe, we need 
an explanation of this change in gas fraction with temperature. 

A number of possible reasons for this trend in gas fraction 
(fgas) could be suggested. These include: 

(1) The change in fgaa is a natural consequence of halo produc- 
tion by stellar mass loss. Any steady state solution would produce 
the trend we see. 

(2) The change in fgas is the product of the surrounding envi- 
ronment. Higher mass systems tend to be in higher mass groups 
and clusters, which have a denser IGM. This prevents the escape of 
gas, or compresses the halo, or allows accretion of gas, causing the 
halo to reach higher densities and gas fractions. 

(3) The lower gas fraction in cooler galaxies is an evolution- 
ary effect, related to the age (or time since last major merger) of 
the galaxy. The mass density of the system is determined by the 
density of the universe at the time at which the system collapsed 
(or underwent its last major merger), so the older systems have the 
highest density and are therefore the hottest. 

(4) Supernova heating, though not the dominant source of en- 
ergy in the halos, does have a second order effect on them. Perhaps 
in lower mass systems, the relatively higher contribution from SNIa 
causes gas to be lost from the potential, lowering the overall frac- 
tion. 

Of these suggested models, the last seems the least likely. Al- 
though it is plausible that supernova heating could have observable 
effects on the galaxy halos, particularly in the least massive mem- 
bers of our sample, it seems unlikely that it could affect the gas 
fraction in the same way at all radii. In the analogous situation of 
gas heating in galaxy groups, the central gas fraction is lowered, 
but only by moving gas to higher radii, raising the gas fraction 
there. There are however some obvious differences which could 
make this comparison invalid, such as the distributed heating and 
mass injection expected from the stars in a galaxy. It is also impor- 
tant to remember that we cannot in most cases observe the galaxy 
halo to high radii. It is possible that gas could be removed to a large 
fraction of the virial radius, where its low density might make it un- 
detectable, and where it could be easily stripped by any interaction 
with a surrounding IGM. However, the fact that we see no sign of 
any trend in fSfu is a strong argument against this hypothesis. 

There are also problems with the idea that such a trend might 
be a natural result of galaxy halos forming from gas of stellar ori- 
gin. Previous work on a larger sample including galaxies found in 
a variety of environments and with a wide range of luminosities 
suggests that galaxy winds are the most likely f ormation method 
for the h alos of early-type galaxies of all sizes iO'Sullivan et alJ 
l2001b) la'). However, these studies also confirmed that large group or 
cluster dominant ellipticals, such as those in our sample, have X- 
ray properties which differ significantly from those of less massive 
galaxies. They could therefore be the most likely group of objects 
to form their gas halos through some other mechanism. Our results 
suggest that the steep trend in the Lx'-Tx relation is caused by the 
drop in fgas and hence Lx- The Lx'-Lb relation also has a steep 
slope compared to that of non-BGG ellipticals iO' Sullivan et alj 
l2001ah . and it seems likely that this is again a product of the trend 
in /gas with temperature. If so, we cannot argue that the trend in 
gas fraction is a product of halo formation by galaxy winds, as we 
would expect all galaxies to show a steep Lx'-Lb relation, as seen 
for this sample. 
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A change in gas fraction caused by the different formation 
histories of the galaxies seems a more feasible model. Assuming 
our galaxies have formed through a number of mergers, we can as- 
sume that there is a scatter in the time since the last major merger. 
While minor mergers can be treated as the accretion of extra matter 
into the galaxy, major mergers will have more profound effects, in- 
cludin g the "resetting" of the po tential well through violent relax- 
ation iSalvador-Sole et anil998l) and the fuelling of bursts of star 
formation. Luminosity weighted ages (based on comparing opti- 
cal absorption line indices to stellar population models) tend to be 
dominated by any young component in the stellar population, so 
that if the galaxy has formed part of its stellar population through 
a merger and starburst, the age measured is likely to be an estimate 
of the time since that merger. In a previous paper fO'SuUi vanetalJ 
we found a relation between the luminosity weighted age of 
a sample of 77 early-type galaxies in a range of environments and 
their normalised X-ray luminosity, Lx/Lb- This indicates that the 
halo luminosity of such galaxies increases as they age. In galaxy 
wind models, this can be explained through the transition from 
rapid outflows to slower winds or inflows, during which the halo 
density and the gas mass in the potential well increases. 

The question then becomes whether we expect the more mas- 
sive galaxies in this sample to have had their last major merger ear- 
lier, and therefore to have higher halo densities. Hierarchical mod- 
els of galaxy formation predict that more massive galaxies have 
longer formation timescales ( Kauffmann 1 996), suggesting that this 
is not the case. However, observational measurements of age for 
sample of galaxies in clusters and the fiel d suggest that mo re mas- 
sive galaxies do indeed form earlier I T homas et alj|2063) . In our 
previous study of X-ray variations with galaxy age we fo und no 
correlation between Lb and age I O' Su llivan et aljflOOld) . but in 
that case our sample included galaxies from a much wider range of 
environments and with more varied luminosities than is the case in 
this paper. In the current sampl e only 16 of our galaxies have mea- 
sured ages in the catalogue of ITerlevich & Forbes 1 200 2), and for 
this subsample we find only a weak ( 1 .4 cr) trend for age to increase 
with L B ■ We are therefore unable to say whether merger history is 
the deciding factor in determining gas fraction, though it does seem 
possible. 

Finally there is the possibility that the change in gas fraction 
we see is related to the environment of the galaxies. The offsets 
we see in the Lx'-Lb and Lb'-Tx relations suggest that the pres- 
ence of a surrounding IGM can increase the X-ray luminosity and 
temperature of the halo. Lx appears to be more strongly affected, 
being offset by a factor of ~9, as compared to a factor of ~2 in 
Tx. A surrounding IGM might cause these offsets by containing or 
stifling galaxy winds, or by through accretion of the IGM on to the 
galaxy halo. If the most massive galaxies were found in the most 
massive groups and clusters we might expect the density of the sur- 
rounding IGM to correlate with galaxy mass. This would mean that 
more massive systems were more strongly affected, and might pro- 
duce the trend in gas fraction with temperature. Of the possible 
processes, containment or stifling would seem the most likely, as 
accretion of IGM gas might be expected to affect gas fraction dif- 
ferently at different radii. One strong argument for this explanation 
is the cor relation between grou p and galaxy X-ray luminosity ob- 
served by Helsdon et al. ('2001). If we assume that this is not caused 
by group cooling flows then it does indeed suggest that galaxy mass 
is correlated with group mass. 



6.4 Environmental Effects and Formation Epoch 

Although most of the relations which we have used to investigate 
our sample of galaxies show no difference between galaxies in dif- 
ferent environments, there are three interesting exceptions. These 
are the Lb'-Tx, Lx'-Lb and jSfit'-Tx relations. The Lb'-Tx and 
Lx'-Lb both show similar behaviour. In each case, we find that the 
best fit relation has a similar slope to that of groups, but is offset to 
higher Lx or Tx at a given Lb - However, when we split the sam- 
ple into galaxies which lie in X-ray bright or X-ray faint groups, 
we find that those in X-ray faint groups do not appear to be off- 
set to the same degree. We have calculated the residuals about our 
best fit lines for both relations, and the residuals we would find if 
the best fit lines for galaxy groups were used. In both cases, the 
mean residual from our best fit line is considerably higher (by 20- 
30%) for the galaxies in X-ray faint groups than for those in X-ray 
bright groups, showing them to be poorly described by the rela- 
tions. The points are almost all offset toward the group relations, so 
these large residuals are a measure of that offset. We can also calcu- 
late the mean residual from the group relations, which will give us 
a measure of the offset of the X-ray faint and bright subsets from 
those relations. The mean residual from the best fit line for groups 
is smaller (by 40-60%) for galaxies in X-ray faint groups. For the 
Lx'-Lb relation, the mean residual of the galaxies in X-ray faint 
groups from the group relation is actually less than half the mean 
residual from the galaxy best fit line, suggesting that although these 
galaxies are not well described by either relation, they are closer to 
the group relation. This supports the idea that the offset is driven by 
the presence of a surrounding IGM. Alternative subdivisions of the 
sample, comparing cooling flow and isothermal galaxies, or group 
dominant galaxies to those in other environments do not show a 
similar segregation, leading us to the conclusion that this difference 
is not a product of group cooling flows. 

There are a number of ways in which a surrounding IGM 
could increase the luminosity and temperature of a galaxy halo. If 
we consider a galaxy in an early wind phase, the rapidly outflowing 
gas might expand until stopped by the pressure of the surrounding 
gas. The IGM could contain the galaxy wind, raising its density and 
preventing the dispersion of high temperature gas to large radii. We 
might therefore expect such an object, even if it has now progressed 
to a later wind phase, to have a higher temperature and luminosity. 
Similarly, if we consider a galaxy in which there is a global inflow 
(within the galaxy halo rather than a group cooling flow initially 
sourced in the surrounding IGM), we might expect gas from the 
IGM to move into the galaxy potential. This gas would add to the 
galaxy luminosity, particularly as the galaxy potential would con- 
centrate it. We might also expect it to have a higher temperature 
than gas produced within the galaxy. In both cases we might also 
see an increased amount of gas in the galaxy at all radii. If the trend 
in gas mass fraction discussed in Section l6. 1.2l is caused by an in- 
crease in halo density for high temperature systems (rather than a 
density decrease in low Tx galaxies), these models might give us 
an explanation of the trend. More massive galaxies would produce 
larger amounts of gas and would be more able to accrete and con- 
centrate a surrounding IGM, so we would expect the trend with 
temperature. 

The other sign of environmental influence on our galaxies is 
shown in Figure|9| If we exclude NGC 1404 (where (ifu is likely 
to be biased by the effects of ram-pressure stripping), we find that 
galaxies in the centres of groups, clusters and cluster sub-clumps 
show a larger degree of scatter in /3/it than those elsewhere. In this 
case the trend seems to be related to the position of the galaxy at 
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the centre of the potential well rather than a surrounding or confin- 
ing IGM. All the non-BGG/BCG galaxies in our sample are found 
in groups or clusters, so their halos are likely to interact with sur- 
rounding gas. A possibility is that the halos of galaxies at the cen- 
tres of larger structures are affected by the gr oup potential, rather 
than the IGM. A recent study of NGC 1399 (Paolillo et al. 2002^ 
shows that in the central regions of that galaxy (the dominant el- 
liptical of the Fornax cluster), the mass profile is dominated by the 
stellar mass of the galaxy. However, the cluster potential does play 
some part at larger radii, and although we have removed the sur- 
face brightness components associated with larger structures, we 
cannot rule out some level of influence. Another possibility is that 
the merger history of these central galaxies has an effect on their 
halo structure. Dynamical friction is thought to cause galaxies in 
groups and clusters to lose orbital energy and fall towards the cen- 
tre of the potential. We might therefore expect to find that group and 
cluster central galaxies undergo many more mergers than galaxies 
elsewhere. The associated disturbance, influx of gas, and star for- 
mation might produce effects which can change the state of the 
galaxy halo, inducing the scatter we observe. 

Although environmental influences seem likely, and provide 
an explanation for the offsets seen in the Lb'-Tx and Lx'-Lb 
relations, this model raises further questions, when compared to 
the other relations. There is a discrepancy between the a:Tx and 
Lb'Tx relations. In the a:Tx relation, temperature appears to be 
a good measure of mass, with a behaviour similar to that seen in 
galaxy clusters. In the Lb'Tx relation we see an offset, suggest- 
ing that while Tx may be a good measure of mass, either galaxies 
have higher temperatures than clusters at a given mass, or that they 
have lower optical luminosities. If we accept the a:Tx relation, we 
must assume that the offset is in Ls , and at first glance this sug- 
gests that galaxies have a higher mass-to-light ratio than clusters. 
This is rather unexpected, as cosmological simulations suggest that 
mass-to-light ratio is correlated with dark halo mass, with galaxy 
clusters havi ng considerably hig her mass per unit L b than individ- 
ual galaxies iSomerville et alli oOl). A similar discrepancy is seen 
between the Lx'Tx and Lx'Lb relations, in which the former has 
an offset to lower Lx (or higher Tx), while the latter has an offset 
to higher Lx (or lower Ls). Once again, if we assume temperature 
behaves as it does in clusters, then the Lx'Tx relation suggests that 
Lx is reduced, relative to larger structures. The offset in the Lx'Lb 
relation would then again imply a large reduction in Ls, compared 
to groups and clusters. It is not clear what is causing these offsets, 
but we can speculate to some extent. 

One possibility to be considered is the difference in formation 
epoch between galaxies, groups and clusters. Although properties 
such as a and X-ray temperature can be considered as measures of 
system mass, they are actually related to the depth of the potential 
well (Tx oc M/R). We expect systems to have densities related to 
the critical density of the universe at the time at which they viri- 
alise (or undergo their last major merger), so systems of a given 
mass which form at earlier epochs should have deeper potentials 
than their more modern counterparts. This suggests that older sys- 
tems will have higher velocity dispersions. X-ray temperatures and 
luminosities than younger ones. A comparison of systems of dif- 
ferent ages should show these differences, and relations involving 
these parameters would be affected them. Relations involving opti- 
cal luminosity might show them particularly well, if the total stellar 
luminosity is related to system mass rather than the depth of the po- 
tential well. 

As we are considering similar systems at different redshifts, 
we take the total mass, Mtot to be constant, but the radius of the 
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system is reduced at higher redshift, R oc (l+z)~^. Therefore, den- 
sity scales as 

Ptot OC OC (1 + z) , (4) 

as expected. We also know that for bremsstrahlung radiation, the 
luminosity can be written 

Lx oc pltR'VT^. (5) 

Since ptotR^ is the mass of gas in the system, we can rewrite this 

as 

Lx oc (1 + zfMtotfgasVT^, (6) 

where fgas is the gas fraction. 

Tx is proportional to M/R, so rearranging Equation^ we can say 
that 

Tx^mI,{1 + z), (7) 
and substituting into Equation|6|will therefore give 

LxOi{l+z)ifgasT%. (8) 

We therefore expect galaxies formed at a redshift Zform = 2 
to have a similar Lx'Tx relation to that of clusters, offset to higher 
luminosities at a given Tx by a factor of 3 2 , and modified by the 
behaviour of fgas- As we in fact see a steeper slope and an off- 
set to lower than expected luminosities, this suggests that fgas is 
both dependent on the Tx of the system (as we have demonstrated 
elsewhere, c.f. Fieure lTTl and in general lower than is the case in 
groups and clusters, countering the expected offset to higher lumi- 
nosities. 

Other relations can also be considered in the same way. Ve- 
locity dispersion scales with mass and radius in a similar way to 
temperature, oc M / R. This means that we expect and Tx 
to scale with redshift by the same factor, so regardless of forma- 
tion epoch the a:Tx relation should remain constant. This agrees 
well with our findings. Optical luminosity is assumed to scale with 
total mass regardless of formation epoch, in which case we can sub- 
stitute it into Equation Q Rearranging the equation shows that we 
would expect a relation of the form Lb oc Tx'^'^, like that of clus- 
ters, but offset to lower luminosities by a factor of 3 2 . This predic- 
tion agrees less well with our results, as we find a slope somewhat 
steeper than that of clusters, with an offset of a factor ~3. We can 
also substitute this predicted relation for Lb'-Tx into Equation|8|to 

find a prediction for the Lx'-Lb relation. This again predicts a re- 

4 

lation unlike that we observe, of the form Lx oc Lb ^ with a large 
offset to high X-ray luminosities. However, this relation contains 
an fgas term which is likely to steepen the relation (as fgas is de- 
pendent on Tx) and partially counter the offset to high Lx- 

6.5 3D models and mass estimates 

Gas fraction is clearly an important quantity for our galaxies. We 
can use the three dimensional models of our galaxies (see Sec- 
tion l5.5> to measure properties such as gas mass and total halo mass 
within a range of radii, in much the same way as we calculate gas 
entropy at a tenth of the virial radius. However it is important to 
emphasise that these models are by necessity inaccurate. We expect 
errors in the virial radii of the galaxies, as we do not know the true 
redshift of formation of each object. We also expect errors associ- 
ated with temperature structure in the galaxy halos, as we assume 
all the halos to be isothermal. The values we find are therefore of 
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more interest as representing trends rather than as exact measures. 
Table|4|hsts mean and median values for our sample. When calcu- 
lating the means, we exclude four galaxies which we find to have 
extreme values of gas fraction. In one case (NGC 4073) we find 
that the gas mass exceeds the mass of the dark halo at all radii, sug- 
gesting that the spatial fit is contaminated by group emission. The 
other three galaxies (NGC 1332, NGC 1549, NGC 4261) all have 
gas fractions of <10~^, considerably lower than most galaxies in 
the sample. Although these are extreme cases, they do demonstrate 
the hazard of accepting measurements from these models at face 
value. The median values we calculate should be unaffected by 
these extreme cases, but the very low median gas fraction and very 
high median mass-to-light ratio within the virial radius make this 
approach somewhat suspect. We have derived the mean gas frac- 
tion, stellar mass fraction and total baryon fraction for the galaxies, 
assum ing a stellar mass-to-light ratio of 5 Mq/Lb iPizzella et al] 
Il997t) . We also assume that the stellar component of the galaxy is 
contained within R^iriai /3. 

The gas masses shown in Table|4]are quite large, as we would 
expect given that most of our sample is made up of groups and clus- 
ter dominant giant ellipticals. iBreeman etalj 11992) find X-ray gas 
masses for their sample of elliptical galaxies in the range ~10*~^^ 
Mq. Our values are comparable with the upper end of this range 
particularly if we take in to account the effect of extrapolation out 
to the virial radius. The mean mass-to-light ratio for th e sample 
may be somewhat high. iLlovd-Davies & PonmanI i2002h find the 
mean mass-to-light ratio of a sample of 20 galaxy groups and clus- 
ters to be 120±20 M ^/Lb, and other comparable estimates in- 
clude 88±33 Mq/Lb (Edge & Stewart 1991) and 100 Mq/Lb for 
the Perseus cluster lEvles et al. 1991) . However, some estimates 
suggest that both elliptical galaxies and gala xy clusters may hav e 
mass-to-light ratios as high as 200 Mq/Lb l lBahcall et ai]|l995l) . 
Taking a conservative view, we assume that the total halo masses 
may be overestimated to some degree, possibly owing to the influ- 
ence of a surrounding group or cluster halo on the galaxies. The 
mean gas, stellar and baryon fractions will all be influenced by any 
overestimate of total mass, so we expect our measured values to be 
lower than is really the case. However, even given such a bias, it is 
clear that the mean gas fraction is considerably lower than is usual 
in galaxy groups and clusters; ~l-2%, as compared to ~20% in 
clusters iLlovd -Davies & Ponman 2002 ; Markevitch & Vikhlinin 
Il997 t). This is exactly the kind of difference we expected to see, 
given the offsets to low luminosity in the Lx'-Tx and Lx'-Lb rela- 
tions. 

The mean star formation efficiency (Mt/[M, -|- Mgaa]) 
is also very high compared to larger structures, which typi- 
cally have e fficiencie s of ~0.2- 0.3 ( Llovd-Davies & Ponman 2002; 
ICirimele et aljil997i ; IXmaud et al. 1992,; David et al. .1990.) . This 
could suggest that galaxies convert much more of their gas into 
stars than is the case in groups or clusters, an unsurprising result 
considering the extent and temperature of a typical cluster halo. 
This result is also consistent with theoretical modelling which sug- 
gests that galaxy sized halos a re likely to c ontain much more cool 
gas that larger systems (Baug h et alll999l) . In ellipticals, cool gas 
is likely to have been formed into stars during the merging pro- 
cess, leading to the high star formation efficiencies we measure. 
However, high star formation efficiencies would also be inferred if 
much of the gas has been removed from the systems via processes 
such as galaxy winds. The mean baryon fraction for our sample is 
quite low compared to measured values for galaxy cluste rs, which 
range from 0.1 1-0.3 iLlovd-Davies & Ponman 2002; Za roubi et alj 
l200lHHradeckv et all200(]l;lDavid et alll995h . The low end of this 



scale is comparable with our value, given the large uncertainties 
we expect in our result. However, for larger values of baryon frac- 
tion in clusters we would have to conclude that up to ^75% of the 
baryons originally in galaxies have been removed, probably being 
blown out into the surrounding IGM. Once again, we emphasize 
that the values listed in Table|4]cannot be taken as precise measure- 
ments. 



7 SUMMARY AND CONCLUSIONS 

We have compiled a sample of 39 massive X-ray luminous early- 
type galaxies for which there are long ROSAT PSPC exposures 
archived. Analysing these data we have carried out detailed spatial 
and spectral fits, and extrapolated from these to approximate three 
dimensional models of the galaxies. The properties measured from 
these fits and models allow us to compare our sample to galaxy 
groups and clusters. Galaxies may be comparable to these larger 
systems because at the simplest level they can be considered as 
dark matter halos containing hot gas, much like the more massive 
structures in which many of them reside. If the dark matter profile 
does not vary significantly with system mass, we might expect the 
gas properties to be similar as well. To compare our galaxies with 
larger galaxy systems we fit a number of relations which are com- 
monly used for groups and clusters, including the Lx:Tx, cr'-Tx, 
Pfit'-Tx, Lb'-Tx and Lx'-Lb relations. We are also able to subdi- 
vide our sample according to environment (galaxies in the centres 
of groups/clusters, galaxies surrounded by a dense IGM) and tem- 
perature structure (galaxies with isothermal halos, central cooling 
or heating). As we have specifically chosen our sample to include 
the most massive galaxies available, it is unsurprising that most of 
them are found to lie in the centres of larger systems and that a 
majority show signs of central cooling. 

We find a number of interesting correlations, and some of the 
relations we fit show similarities to those found for groups and clus- 
ters. These include; 

(1) The a:Tx relation is identical (within errors) to that found 
for galaxy clusters, suggesting that both gas temperature and cen- 
tral velocity dispersion are good measures of the depth of the po- 
tential well. In clusters this correlation arises because gas is heated 
through shocks as it falls into the potential and this may also be 
the case for our galaxies. However, it also seems possible that tem- 
perature could be related to the depth of the potential if the gas is 
mainly of stellar origin, at least in large galaxies such as those in 
our sample. 

(2) The Lx :Tx relation has a very similar slope to that found for 
galaxy groups. Although this could be caused if group scale cooling 
flows are the main source of the halos of our sample of galaxies, 
the galaxies in which we find evidence of cooling do not behave 
differently from those which show no signs of cooling, in this or 
any of our other relations. This suggests that cooling flows are not 
the predominant source of emission in our sample of galaxies. 

(3) The Lx'-Lb and Lb'-Tx relations also have similar slopes to 
those found for galaxy groups. 

(4) The X-ray surface brightness profiles of the galaxies in the 
sample have a mean [3 fa parameter of ~0.5, similar to the mean 
optical surface brightness slope expected for elliptical galaxies. 
This suggests that there is equipartition of energy between the gas 
and stellar components of these systems, a conclusion we also draw 
from the slope of the a'-Tx relation. However, there is a variation in 
the normalisation of the profiles with temperature, indicative of a 
variation in the overall gas fractions of the systems. This variation is 



© 0000 RAS, MNRAS OOQ.ITlll^ 



X-ray Scaling Properties of Early-type Galaxies 25 



Mean Median 





Rvirial 


Rvirial / 3 


f^virial 


f^virial / 3 


Total mass (Mq ) 


1.31x1013 


4.25x1012 


1.21 xlOi^* 


3.60x1012 


Gas mass (Mq) 


2.39x10" 


4.60x101" 


6.26x10"' 


1.89x101° 


Stellar mass (Mq) 


4.75x1011 


4.75x1011 


3.53x1011 


3.53x1011 


Gas fraction 


0.019 


0.012 


0.006 


0.004 


Stellar fraction 


0.050 


0.16 


0.038 


0.113 


Baryon fraction 


0.069 


0.17 


0.049 


0.124 


Star Formation Efficiency 


0.77 


0.92 


0.84 


0.96 


Mass-to-Light ratio (MqILb) 


166.7 


54.8 


281.1 


38.62 



Table 4. Mean and median masses and other derived quantities for our sample of galaxies. 



likely to be the cause of the steep slope of the Lx'-Tx relation, with 
lower temperature systems having lower gas fractions and therefore 
lower X-ray luminosities. 

We also find a number of differences between the relations 
found for galaxies and those known to exist in groups and clusters, 
such as: 

(1) The Lx'-Tx, Lk'-Lb and Lb'-Tx relations, which have sim- 
ilar slopes to the relations found for galaxy groups, are offset from 
those relations. These offsets could be at least partially explained 
if the halos of the galaxies were in fact group cooling flows, but as 
mentioned above this appears to be unlikely. 

(2) There is no trend in mean gas entropy with temperature, 
which we might expect to find if galaxy halos were formed in the 
same way as those of galaxy groups and clusters. 

(3) There is also no trend in surface brightness slope with tem- 
perature. In groups this trend is caused by the heating and move- 
ment of gas to higher radii; in galaxies our results suggest that this 
does not occur. Instead we find a trend in overall gas fraction with 
temperature. Could this be a sign that heating expels gas from the 
system, or that the halo builds up to higher densities and tempera- 
tures over time? 

(4) In the Lb'-Tx and Lx'-Lb relations we find some evidence 
that galaxies in groups whose X-ray halos are poor do not fol- 
low the same relations as those galaxies which are surrounded by 
a rich IGM. This may be evidence of environmental influence on 
our galaxies, though it does not appear that this effect is caused by 
cooling flows. 

Taken as a whole, our results seem to suggest that although 
galaxy X-ray halos obey relations similar to those found for groups 
and clusters, they are not formed through the same process. It seems 
likely that a significant portion of the gas in the halos of early-type 
galaxies is produced by stellar mass loss within the galaxy, rather 
than infall and shock heating of primordial material. The similari- 
ties in the behaviour of the halos across the range of system sizes 
would in this case be a product of the dark matter potential, which 
is believed to follow a profile of similar form, regardless of sys- 
tem mass. Factors such as the formation and merger history of the 
galaxy, gas cooling and the surrounding environment are all likely 
to influence the development of the halo, and further exploration 
of the properties of early-type galaxies using higher quality X-ray 
data may allow some of the outstanding issues raised by this study 
to be resolved. It is particularly important to remember that a large 
proportion of the galaxies in this sample lie at or near the centres 
of galaxy groups. The advent of Chandra and XMM-Newton makes 
an improved study in this area possible, once a sufficient number 
of observations of X-ray bright early-type galaxies in a range of 
environments become available. 
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